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ABSTRACT 





This thesis investigates the potential use of wurtzite 








Indium Gallium Nitride as photovoltaic material. Silvaco 





Atlas was used to simulate a quad-junction solar cell. Each 








of the junctions was made up of Indium Gallium Nitride. The 





band gap of each junction was dependent on the composition 











percentage of Indium Nitride and Gallium Nitride within 





Indium Gallium Nitride. The findings of this research show 








that Indium Gallium Nitride is a promising semiconductor for 





solar cell use. 
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EXECUTIVE SUMMARY 











One of the primary goals of solar cell design is to 
improve efficiency. Indium Gallium Nitride is a material 
that has undergon xtensiv research since 2002 as a 
potential photovoltaic material. By varying the composition 














of Indium Nitride and Gallium Nitride within Indium Gallium 








Nitride, the band gap of this semiconductor material can be 








changed. The band gap range of Indium Gallium Nitride 
matches closely the visible solar spectrum frequencies. 


Hence, a high-efficiency solar cell can be potentially 





developed by having several Indium Gallium Nitride 


junctions. 


The use of Silvaco Atlas as a simulation tool can aid 








in identifying whether Indium Gallium Nitride can be used 





for photovoltaics. Previous research at the Naval 
Postgraduate School has demonstrated the viability of 


Silvaco Atlas as a solar cell modeling mechanism. 


Finding high-efficiency solar cell models is of great 


interest in space applications. By increasing the efficiency 





of photovoltaics, the number of solar panels is decreased. 





Therefore, the overall weight that needs to be launched into 





space is reduced. A cost-reduction can be achieved over th 


life of the space power system. 


Research at lLawrenc Berkeley National Laboratory 











(LBNL) has been ongoing since 2002 in trying to develop 





high-efficiency solar cells. One of the materials LBNL has 








been investigating is Indium Gallium Nitride. In addition to 














the band gap range, Indium Gallium Nitride also has other 








characteristics that are beneficial for space applications. 
XV 





The radiation level that Indium Gallium Nitride is able to 








withstand is approximately two orders of magnitude greater 


than current multijunction solar cell materials. Therefore, 











it is of great interest to continue to research Indium 








Gallium Nitride and its possibilities as photovoltaic 


material. 


The goal of this research was to investigate the 











efficiency of Indium Gallium Nitride solar cells using the 








Silvaco Atlas TCAD simulation software. No other TCAD 





simulations have been performed on InGaN solar cells. The 
methodology of this research progressed from single to 
multi-junction cells. The simulation results predicted 


efficiencies as high as 41% for a four-junction solar cell. 





Further refinements of the simulation model are. still 











possible. Actual production of a single-junction solar cell 








is the next step required to eventually manufacture high- 











efficiency, multi-junction InGaN solar cells. 


XV1 


I. INTRODUCTION 


A. BACKGROUND 


One of the most important problems to solve in space 


solar cell design is efficiency. As efficiency increases, 





the required number of cells decreases to fulfill electrical 





power requirements. 


Solar cell efficiencies have improved over time by 
increasing the number of junctions. Each junction is capable 


of extracting energy from a portion of the solar spectrum. 


A new path to improving efficiency is to use new 





photovoltaic materials such as wurtzite Indium Gallium 














Nitride. Indium Gallium Nitride possesses a band gap range 
that can extract energy from a large portion of the visible 


solar spectrum. 


B. OBJECTIVE 





The objective of this thesis is to simulate Indium 





Gallium Nitride solar cells to predict their efficiencies. 





After obtaining the efficiencies, a comparison is made with 

















current solar cell efficiencies. Silvaco Atlas TCAD 
simulation software is used to investigate these 
efficiencies. 
Cc. RELATED WORK 

Previous Silvaco Atlas solar cell simulations have been 
performed by Naval Postgraduate School researchers. 
Michalopoulos [1] investigated the feasibility of designing 
solar cells using Silvaco Atlas. To demonstrate the use of 








1 


Silvaco Atlas, Michalopoulos simulated single-junction solar 








cells with Gallium Arsenide, dual-junction solar cells with 














Indium Gallium Phosphide and Gallium Arsenide, and triple- 











junction cells with Indium Gallium Phosphide, Gallium 
Arsenide, and Germanium. The highest efficiency obtained 


with the triple-junction was 29.5%. This result matched the 





29.3% efficiency obtained with actual triple-junction solar 
cells in production. Bates [2] provided excellent background 
on the use of Silvaco Atlas. Green [3] and Canfield [4] 


continued to work on Silvaco Atlas solar cell design. Matlab 





functions and Silvaco Atlas code were modified from [1]-[4] 


to perform simulations in this research. 





A previous Indium Gallium Nitride simulation has been 








performed Bon using fundamental semiconductor physics 
formulas. The results of this research are compared with 


those OF |[:5:)). 


D. ORGANIZATION 


1. Purpose of Solar Cells 


The purpose of solar cells is to convert light energy 
into electrical energy. Light is made up of photons. Photons 


carry energy that is dependent on the color, or wavelength 





of light. Electrical energy is generated when photons excite 





electrons from the valence band into the conduction band in 











semiconductor materials. Chapter covers fundamental 





concepts of solar cells and semiconductor materials. 
2. Simulation Software 
Silvaco Atlas was used as the simulation software in 


this thesis. Silvaco offers a suite of software programs 
2 


that can act as a Virtual Wafer Fabrication (VWF) tool. The 


Naval Postgraduate School has had multipl researchers 








publish theses on the use of Silvaco Atlas for the purpose 

















of solar cell modeling. Chapter covers the basics of 


Silvaco Atlas. 


3. Indium Gallium Nitride 








Wurtzite Indium Gallium Nitride (InGaN) is a 


semiconductor that has the potential to produce high- 





efficiency solar cells. Dr. Wladek Walukiewicz, of Lawrence 
Berkeley National Laboratory (Solar Energy Materials 


Research Group), has been conducting extensive InGaN 





research for the purposes of producing photovoltaic 








material. Chapter IV covers the characteristics of InGaN, 


its potential, its material science research status, and its 











current limitations. 
4. Simulation 
Based on the optical characteristics of InGaN, 


simulations were performed using Silvaco Atlas. Chapter V 





covers the findings of the simulations. Chapter VI covers 








the conclusions and recommendations. Finally, the appendices 





include the code used in Silvaco Atlas as well as Matlab. 
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II. SOLAR CELL AND SEMICONDUCTOR BASICS 





This thesis examines a novel process to simulate Indium 





Gallium Nitride solar cells. Before delving into the 


Simulation, this chapter covers the fundamental principles 





of solar cells and semiconductors. The properties of the 


semiconductor material determine the characteristics of the 





photovoltaic device. 


A. SEMICONDUCTOR FUNDAMENTALS 


1. Classification of Materials 


Materials can be categorized according to their 
electrical properties as conductors, insulators or 
semiconductors. The conductivity o is a key parameter in 


identifying the type of material. 


Figure 1 presents a sample of materials based on 
conductivity. The semiconductors fall between the insulators 


and the conductors. 


Resistivity p (Q — cm) 
108 106 yo 10! 10! 108 108 104 10” 1 102 104 10% 104 


Germanium (Ge) Silver 
© Glass — 
Nickel oxide Silicon (Si) Copper 
(pure) 


i sallium : i 3 Aluminum 
Diamond Gallium arsenide (GaAs) ni 


(pure) 


Gallium phosphide (GaP) Plotinum 


© Sulfur 


Fused Cadmium sulfide (CdS) Bismuth 


quartz 





108 106 jo jo? Jo 10% 1o* 104% 107% 1 10? —:104 10° ~— 108 
Conductivity o (S/cm) 


| Insulator ——_+|_— Semiconductor ——+— Conductor —| 


Figure 1. Materials classified by conductivity (From [6]). 





5 


Semiconductors can be found in elemental or compound 
form. Silicon (Si) and Germanium (Ge) are examples of 


elemental semiconductors. Both of these semiconductors 





belong to group IV of the periodic table. 








Figure 2 shows an abbreviated periodic table. In 





addition to the group IV semiconductors, compounds can be 














made with elements from groups and V, respectively. 




















Examples of -V semiconductors include Aluminum Phosphide 








(ALP), Gallium Nitride (GaN), Indium Phosphide (InP), 








Gallium Arsenide (GaAs), among others. It is also possible 














to make semiconductor compounds from groups -VI, such as 





Zinc Oxide (ZnO), Cadmium Telluride (CdTe), Mercury Sulfide 


























(HgS), among others. The subject of this thesis is the 
ternary alloy wurtzite Indium Gallium Nitride (InGaN). 
Indium and Gallium are group elements, while Nitrogen is 


a group V element. 





Figure 2. Partial periodic table (After [7]). 





2. Atomic Structure 


Since Silicon is the most commonly used semiconductor 
in solar cells today, a brief analysis of this element is 


presented. 


Silicon has 14 protons and 14 neutrons in its nucleus. 


The 14 electrons are distributed in three shells. 


Figure 3 shows the arrangement of electrons in a 


Silicon atom. There are two electrons in the first shell, 





eight electrons in the second shell and four electrons in 





the outer shell. When Silicon atoms are together, the atoms 





from the outer shells form covalent bonds. Hence, a Silicon 


atom forms bonds with four other Silicon atoms. 





Figure 3. Silicon atomic structure (From [8]). 


As shown in Figure 4, the covalent bonds formed by the 








silicon atoms ar represented by the elliptical dotted 





lines. In the absence of an electron, the covalent bond 





ceases to exist. A hole takes the place of the electron. 
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Figure 4. Silicon atom covalent bonds (From [1l, p. 11]). 


Energy bands are a fundamental concept of semiconductor 
physics. These energy bands are the valence band, the 
conduction band, and the forbidden gap or band gap [9, p. 


27]. When an electron is in the valence band, the covalent 





bond exists. In order for the electron to move from the 
valence band into the conduction band, energy is required to 


excite the electron. The band gap energy is the minimum 








energy required for the electron to make the move from the 


valence band into the conduction band. 


Figure 5 shows the band gap diagrams for conductors, 
insulators, and semiconductors. This figure also reinforces 


the differences among the three types of materials according 








to their electric properties. In the case of conductors, the 
band gap is small or non-existent. By contrast, insulators 
have wide band gaps. Therefore, it takes much more energy 


for the insulator to have electrons in the conduction band. 










Bandgap ! 1 Bandgap Ez 


Valence band 


Valence band 





Valence band 
Conductors Semiconductors Insulators 
Figure 5. Band gap diagrams (From [1l, p. 8]). 


Semiconductors have band gaps that are dependent on the 


material. Because band gap is also dependent on temperature, 





it should be noted that all quoted band gaps in this thesis 


are at room temperature (300 K.) 


3. Electrons and Holes 


A semiconductor at absolute zero temperature is unable 





to conduct heat or electricity [10, p. 44]. All of the 
semiconductor’s electrons are bonded. The electrons acquire 
kinetic energy as the temperature is increased. Some of the 


electrons are freed and move into the conduction band. These 





electrons are able to conduct charge or energy. The areas 
left by these electrons in the valence band are called 


holes. As the temperature of the semiconductor increases, 








the number of free electrons and holes increases as well. 











Therefore, the conductivity of the semiconductor is directly 


proportional to temperature increases. 


The conductivity of the semiconductor can be altered by 
exposing it to light or by doping it. Photoconductivity 


consists of exposing the semiconductor with photon energy 





larger than the semiconductor’s band gap. Doping consists of 


adding impurities to the semiconductor material. 





Figure 6 shows n-type and p-type doping. The atom of 
the semiconductor is represented by the blue circle and a 
+4. For example, a silicon atom has four electrons in its 
outer shell. When silicon atoms are n-doped, atoms from 
group V of the periodic table are added. Each added group V 


atom provides an extra electron to donate. When silicon 














atoms are p-doped, the atoms from group of the periodic 














table are added. Each group atom has one less electron 











than silicon. Therefore, a hole is added. Doping increases 


the number of carriers in the semiconductor material. 
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Figure 6. Doping: n-type and p-type (From [1, p. 13]). 


An intrinsic semiconductor is a pure semiconductor with 
a negligible amount of impurity atoms [9, p.34]. By 
definition, the number of electrons and holes in a 


semiconductor are represented by n and p. See Table 1. 
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_ number of electrons 


3 
cm 


_ number of holes 


cm? 








Table 1. Definitions of n and p. 





In the case of an intrinsic semiconductor, the 





following case occurs: 


n= p=ni 


For example purposes, ni is given for GaAs, Si and Ge 








and room temperature [9, p.34]. Data for InN and GaN are 


from [11]. See Table 2. 
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Table 2. n;, for five semiconductors. 
When solar cells utilize semiconductor materials from 
Table 2, the current level is highest for Ge and lowest for 





GaAs. The current levels correspond in rank to the n; levels 








provided in Table 2. No data exists for InN or GaN current 


levels as photovoltaic materials. 





4. Direct and Indirect Band Gaps 


Since the band gap is the minimum energy required to 
move an electron from the valence band into the conduction 
band, it is necessary to distinguish between direct and 


indirect band gaps. 


Figure 7 shows the concept of direct and indirect band 
gaps. The blue portion represents the valence band. The tan 


portion represents the conduction band. 








Figure 7. Direct and indirect band gaps (After [1l, p. 24]). 





When the valence band and the conduction band coincide 


in wave vector k, the semiconductor has a direct band gap. 





When the valence band and the conduction band have different 
wave vector k, the semiconductor has an indirect band gap. 


The k vector represents a difference in momentum. Photons 





have negligible momentum. In order to excite an electron 
from the valence band to the conduction band in an indirect 
band gap semiconductor, in addition to a photon, a phonon is 


required. The phonon is a lattice vibration. The phonon 
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transfers its momentum to the electron at the time the 





photon is absorbed. Therefore, a direct band gap 
semiconductor is generally better for optoelectronics. 


Silicon is an example of an indirect band gap semiconductor. 








Gallium Arsenide and wurtzite Indium Gallium Nitride are 


examples of direct band gap semiconductors. 


53 Fermi Level 


The Fermi function f(E) specifies how many of the 
existing states at energy E are filled with an electron [9, 


p.42]. The Fermi function is a probability distribution 





function defined as follows: 


f(E)= ta 


Ite ‘FT 





Where Erp is the Fermi level, k is Boltzmann’s constant, 





and T is the temperature in Kelvin. 


From the Fermi function, it can be determined that when 





E=Er, then f (E)=f (Erp) =0.5. 


Figure 8 shows the Fermi level for an intrinsic 


semiconductor. Figure 9 shows the Fermi level for n-type 





semiconductor. Figure 10 shows the Fermi level for p-type 








semiconductor. From these Figures, it can be deduced that 


the n-type material has a larger electron carrier 








distribution and th p-typ material has a larger hole 


carrier distribution. 
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Figure 8. Fermi level: intrinsic case (After [9, p. 42]). 
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Figure 9. Fermi level: n-type case (After [9, p. 42]). 
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Figure 10. Fermi level: p-type case (After [9, p. 42]). 
B. SOLAR CELL FUNDAMENTALS 


After covering the basics of semiconductors, the 


logical step is to continue with solar cell fundamentals. 


1. History of Solar Cells 


Solar cell research started when Edmund Bequerel 





discovered the photovoltaic effect in 1839. He recorded that 








an electric current was produced when light was applied to a 
silver coated platinum electrode immersed in electrolyte. 


The next significant step was performed by William Adams and 





Richard Day in 1876. They discovered that a photocurrent 





appeared when selenium was contacted by two heated platinum 
contacts. However, current was produced spontaneously by the 


action of light. Continuing to build on these efforts, 





Charles Fritts developed the first large area solar cell in 
1894. He pressed a layer of selenium between gold and 
another metal. Progress continued as the theory of metal- 
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semiconductor barrier layers was established by Walter 





Schottky and other. In the 1950s, silicon was used for solid 
state electronics. Silicon p-n junctions were used to 
improve on the performance of the Schottky barrier. These 


silicon junctions had better rectifying action and 





photovoltaic behavior. In 1954, Chapin, Fuller and Pearson 
developed the first silicon solar cell, with a reported 
efficiency of 6%. However, the cost per Watt associated with 


these solar cells made them prohibitively xpensiv for 








terrestrial use. However, locations where power generation 








was not feasible (i.e., space) were suitable for solar 


cells. Satellites wer th first clear application for 








silicon solar cells. Since that time, solar cells have 
progressed steadily both in terms of efficiency as well as 


materials used their production [10, p. 2]. 


Table 3 shows a brief list of events in the history of 


photovoltaics from 1939 until 2002. 





2. The Photovoltaic Effect 


The photovoltaic effect is the process by which a solar 








cell converts the energy from light into electrical energy. 
Light is made up of photons. The energy of these photons 
depends on the color (wavelength) of light. The material 
that makes up the solar cell determines the photovoltaic 


properties when light is applied [10, p. 1]. 


When light is absorbed by matter, such as metal, 
photons provide the energy for electrons to move to higher 


energy states within the material. However, the excited 











electrons return to their original energy state. In 


semiconductor materials, there is a built-in asymmetry (band 


IS 


gap). This allows the electrons to be transferred to an 





external circuit before they can return to their original 





energy state. The energy of the excited electrons creates a 
potential difference. This electromotive force directs the 


electrons through a load in the external circuit to perform 





electrical work. 


e 1839 Becquerel (FR) discovered photogalvanic effect in liquid electrolytes 

e 1873 Smith (UK) discovered photoconductivity of solid Se 

e 1877 Adams and Day (UK) discover photogeneration of current in Se tubes; the first observation of 

PV effect in solids 

1883 Fritts (US) makes first large area solar cell using Se film 

1954 First 6% efficient solar cells reported: Si (Bell Lab, USA) and CugS/CdS (Air Force, USA) 

1955 Hoffman Electronics (USA) offers 2% efficient Si PV cells at $1500/W 

1958 NASA Vanguard satellite with Si backup solar array 

1959 Hoffman Electronics (USA) offers 10% efficient Si PV cells 

1963 Sharp Corp (JP) produces first commercial Si modules 

1966 NASA Orbiting Astronomical Observatory launched with | KW array 

1970 First GaAs heterostructure solar cells by Alferov, Andreev et al. in the USSR 

1972 First PV conference to include a session on terrestrial applications (IEEE) 

1973 A big year in photovoltaics: Worldwide oil crisis spurs many nations to consider renewable 

energy including photovoltaics; Cherry Hill Conference in USA (established photovoltaics’ potential 

and legitimacy for government research funding); World's first solar powered residence (University of 

Delaware, USA) built with Cu)S (not c-Si!) solar modules 

e 1974 Project Sunshine initiated in Japan to foster growth of PV industry and applications; Tyeo (USA) 

grows 2.5 cm wide Si ribbon for photovoltaics, first alternative to Si wafers 

1975 First book dedicated to PV science and technology by Hovel (USA) 

1980 First thin-film solar cell > 10% using CuyS/CdS (USA) 

1981 350 KW Concentrator array installed in Saudi Arabia 

1982 First | MW utility scale PV power plant (CA, USA) with Arco Si modules on 2-axis trackers 

1984 6 MW array installed in Carrisa Plains CA, USA [35] 

1985 A big year for high-efficiency Si solar cells: Si solar cell >20% under standard sunlight 

(UNSW, Australia) [36] and >25% under 200X concentration (Stanford Univ. USA) [37] 

1986 First commercial thin-film power module, the a-Si G4000 from Arco Solar (USA) 

e 1987 Fourteen solar powered cars complete the 3200 km World Solar Challenge race (Australia) with 
the winner averaging 70 kph 

e 1994 GalnP/GaAs 2-terminal concentrator multijunction > 30% (NREL, USA) [38] 

e 1995 “1000 roofs” German demonstration project to install photovoltaics on houses, which triggered 
the present favorable PV legislation in Germany, Japan and other countries 

e 1996 Photoelectrochemical “dye-sensitized” solid/liquid cell achieves 11% (EPFL, Switzerland) [39] 

e 1997 Worldwide PV production reaches 100 MW per year 

e 1998 Cu(InGa)Se, thin-film solar cell reaches 19% efficiency (NREL, US) [40] comparable with 
multicrystalline Si. First concentrating array for space launched on Deep Space | by US (5 kW using 
high efficiency GalnP/GaAs/Ge triple junction cells) 

e 1999 Cumulative worldwide installed photovoltaics reaches L000 MW 

e@ 2000 Olympics in Australia highlight wide range of PV applications, and the awarding of the first 
Bachelor of Engineering degrees in Photovoltaics and Solar Engineering (UNSW, Australia) 

e@ 2002 Cumulative worldwide installed photovoltaics reaches 2000 MW. It took 25 years to reach the 
first 1000 MW and only 3 years to double it; production of crystalline Si cells exceeds 100 MW per 
year at Sharp Corp. Japan). BP Solar ceases R&D and production of a-Si and CdTe thin-film modules 
in USA ending >20 years of effort 


Table 3. Notable events in the history of photovoltaics (From 
[12]). 
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a. The Electromagnetic Spectrum 


Light is electromagnetic radiation. The frequency 





of light determines its color. Figure 11 shows the visible 
part of the electromagnetic spectrum. Visible wavelengths 


range from 390 nm (violet) to 780 nm (red). 





VISIBLE 
Wavelength (m) 
w? ws 107 6 105 10" 
Ultraviolet 
| 
Frequency (Hz) 
103 102 10" 
Photon energy (eV) 
SOLAR 
RADIATION 
Figure 11. The electromagnetic spectrum (From [13]). 


Table 4 shows the approximate wavelength range of 


visible colors. 





























Color Wavelength (nm) 

Red 780-622 

Orange 622-597 

Yellow 597-577 

Green 577-492 

Blue 492-455 

Violet 455-390 

Table 4. Approximate wavelength of various colors in vacuum 

(After [14]). 
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The sun emits light from ultraviolet, visible, and 
infrared wavelengths in the electromagnetic spectrum. Solar 


irradiance has the largest magnitude at visible wavelengths, 





peaking in the blue-green [10, p. 17]. 


b. Band Gap 


The band gap of the semiconductor material 


determines how the solar cell reacts to light. Table 5 shows 





a small sample of semiconductor band gaps. Chapter IV covers 














the band gaps of Indium Nitride, Gallium Nitride, and Indium 


Gallium Nitride. 
































Material Band gap (eV) at 300 K 
Si Lele? 
Ge 0.66 
GaAs 1.42 
InP 1.34 
Table 5. Common semiconductor band gaps (After [9, p. 31]). 


The band gap of the semiconductor material 





determines the wavelength of light that meet the 








requirements to generat lectrical energy. The conversion 








formula between band gap and wavelength is: 


he 
Mee) = eV) 

1.24 
Mum) == (eV) 


Where A is the wavelength in micrometers, h is 


Planck’s constant, c is the speed of light in vacuum, and Eg 
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is the band gap in eV. One eV iS approximately equal to 





1.6x10779 J of energy. In the case of Gallium Arsenide, the 





wavelength that corresponds to 1.42 eV is 0.873 mum. 


Figure 12 helps visualize the concept of light 
absorption. When light has energy greater than 1.1 eV, the 


Silicon solar cell generates electricity. Light with less 





than 1.1 eV of energy is unused. Similarly, light with 





energy greater than 1.43 eV excites the outer. shell 
electrons of the gallium arsenide solar cell. And finally, 
light with energy greater than 1.7 eV is useful for aluminum 


gallium arsenide photovoltaic material. 








>ttev 





<I.dev < 143 eV <1.7ev 
Figure 12. Effect of light energy on different band gaps 
(From [15]). 
c. Solar Cell Junctions 


The discussion in the previous section discussed 
the effect of light energy on different band gaps. Treated 


individually, each of the photovoltaic materials from Figure 





12 would act as a single junction solar cell. 


However, to increase the efficiency of the solar 





cell, multiple junctions can be created. For example, in 
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Figure 12, the top junction would be made up of Aluminum 
Gallium Arsenide. This junction would absorb light energy 


greater than 1.7 eV. Any unused photons would be filtered 





through to the next junction. The gallium arsenide junction 
would then absorb the photons with energy greater than 1.4 
eV. The remaining photons would be absorbed by the silicon 


junction. 


Although the above paragraph described the basics 
of a multijunction solar cell, such device may not produce 


the desired results due to lattice mismatch. The next 





section covers the basics of lattice matching. 


d. Lattice Matching 





Semiconductors are three-dimensional in their cell 





structure. The simple cubic structure serves to illustrate 





the concept of lattice and is presented in Figure 13. 


a 


Figure 13. Simple cubic lattice structure (From [16]). 





Figure 13 shows that each side of the cube is 
represented by the letter “a”. The separation “a” is called 
the lattice constant. Each semiconductor material has a 
lattice constant. Therefore, when creating multijunction 


solar cells, the lattices must be matched. 
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Figure 14 shows the lattice constants for several 


semiconductors. An example given by P. Michalopoulos [1, p. 








87] shows how to lattice match Indium Gallium Phosphide to 


Gallium Arsenide. 
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Figure 14. Lattice constants (From [17]). 





The alloy Indium Gallium Phosphide is composed of 





x parts of Gallium phosphide and 1-x parts of Indium 





Phosphide. Therefore, Indium Gallium Phosphide is 





represented as Inj-,Ga,P. 


GaAs has a lattice constant a=5.65A, GaP has 





a=5.45A and InP has a=5.87A. The goal is to create InGaP 





with a lattice constant that matches that of GaAs. The 


formula to find x is given as: 
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a —-aA 
oes — “GaAs InP 
BGaas =AGap Xt App U-x) x= 


AGap ~ np 


With x=0.52, Ino.agGao.52P has o=5.65A. A rough 








approximation of the resulting band gap of Ino.4gGao.52P is 


given as follows: 
InGaP _ GaP InP 
EnoaP — EC”? + (1- EB 


The equation yields a band gap of 1.9 eV for 








Ino.4gGao.52P. Therefore, a dual-junction solar cell of InGaP 


at 1.9 eV and GaAs at 1.4 eV can be built. 





Figure 15 shows the lattice constant for Indium 


Nitride and Gallium Nitride. 
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Figure 15. Lattice constant for InN and GaN (From [18]). 


e. AMO Spectrum 


The location of the solar cell affects the input 


solar radiation spectrum. A solar cell on Mars receives a 
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different (smaller) spectrum than a solar cell on ia 
satellite that orbits Earth. The energy received outside 
Earth’s atmosphere is approximately 1365 W/m’. This spectrum 


is called Air Mass Zero or AMO. Terrestrial solar cells have 





to deal with the attenuation of the solar spectrum due to 


EBarth’s atmosphere. This solar spectrum is called AM1.5. For 








the purposes of this thesis, AMO is used during simulations. 





From data obtained from the National Renewable Energy 
Laboratory (NREL) [19], the AMO spectrum was plotted using a 
Matlab script. 


Figures 16 and 17 show the AMO spectrum with 
respect to wavelength and energy, respectively. From Figure 
17, it can be seen that semiconductor materials with band 


gaps of less than 4 eV are able to extract most of the solar 








spectrum. 
AMO 
0.25 
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Wavelength (um) 
Figure 16. AMO spectrum (Wavelength vs Irradiance) (After 


[19]).. 
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Figure 17. 
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AMO spectrum (Energy vs Irradiance) (After [19]). 
£. Current-Voltage Curves 
A typical solar cell current-voltage (IV) curve is 





presented in Figure 18. 


Figure 18. 
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[20]). 





From Figure 18, ther ar several points of 





interest. The short circuit current (Isc) occurs when the 


voltage is zero. This is the highest absolute value current. 





The open circuit voltage (Voc) occurs when the current is 


zero. This is the highest voltage. The dimensions of the 





larger rectangle in Figure 18 are determined by Voc and Isc. 





Since power (P) is determined by the product of current 





times voltage, the maximum power point occurs at (Vyp,Iyp) . 


The calculations for solar cell efficiency are as 





follows: 


Pmax=Imp V mp 
= Pimax = Imp Vmp 
Isc Voc Isc Voc 
Pimax = Imp V mp 
Pin Pin 








7) = 


Where Pmax is the maximum power point, FF is the 


fill factor, and 9 is the efficiency. The fill factor 





measures the “squareness” of the IV curve. 





g. Electrical Output 





A solar cell iS a p-n junction photodiode. In 








order to obtain the IV characteristic of the solar cell, the 





dark current needs to be subtracted from the photogenerated 





current. 
l=t.-Ib 


The dark current is the current through the solar 


cell when bias is applied in the dark [10, p. 30]. 





Graphically, the IV characteristic is obtained in 


Figure 19. 
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V Vv 
Figure 19. Solar cell IV characteristic (From [21]). 
Cc. CHAPTER CONCLUSIONS 
This chapter provided basic information on 


semiconductors and solar cells. The foundation has been 
established from the physics stand point. The next step is 


to cover an introduction to the simulation software. 
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IIIT. SILVACO ATLAS SIMULATION SOFTWARE 


This thesis uses Silvaco Atlas to perform solar cell 





modeling. Silvaco International produces a suite of software 


programs that together become a Virtual Wafer Fabrication 





tool. This chapter introduces Silvaco Atlas and some of its 


features. 


A. VIRTUAL WAFER FAB 





Silvaco International provides several software tools 





to perform process and device simulation. From Figure 20, it 


can be seen that Silvaco offers powerful simulation 





software. 


ATHENA 


MERCURY 
FastFET Simulation 
Software 


SSuprem3 ATLAS Pagsica based UTMOST 


1D/2D Process 
Simulation Software 


1D Process Device Simulation 
Simulation Software Software 


SPICE Modeling 
Software 





Transport Parameter Parasitic Extraction 
Generator Software 


bf bt JT Yt UT Ut Sb 


Virtual Wafer Fab Interactive Environment 
Utility Layout Editor 


Figure 20. Silvaco’s Virtual Wafer Fabrication Environment 
(From [22]). 


DeckBuild DevEdit TonyPlot 1D/2D TonyPlot 3D C-Interpreter 
Run-Time Structure & Mesh [| Interactive Visualization (J Interactive Visualization J User Defined 
i ili Utili 


Environment Editor Utility ility Models 





Platforms: Solaris, HPUX, Linux and Windows 








In this thesis, Silvaco Atlas was extensively used. The 
DeckBuild run-time environment received the input files. 


Within the input files, Silvaco Atlas was called to execute 





the code. And finally, TonyPlot was used to view the output 
of the simulation. Additionally, output log files were 


produced. The data extracted from the log files could then 
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be displayed using Microsoft Excel or Matlab scripts. Figure 


21 shows the inputs and outputs for Silvaco Atlas. 


C DevEdit 
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Figure 21. Atlas inputs and outputs (From [23, p. 2-2]). 


B. SILVACO ATLAS 


Atlas is a software program used to simulate two and 


three-dimensional semiconductor devices. The physical models 








include in Atlas are presented in Table 6. 
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e DC, AC small-signal, and full time-dependency. 
e Drift-diffusion transport models. 


e Energy balance and Hydrodynamic transport models. 


e Lattice heating and heatsinks. 
e Graded and abrupt heterojunctions. 


e Optoelectronic interactions with general ray tracing. 


e Amorphous and polycrystalline materials. 

e General circuit environments. 

e Stimulated emission and radiation 

e Fermi-Dirac and Boltzmann statistics. 

e Advanced mobility models. 

e Heavy doping effects. 

e Full acceptor and donor trap dynamics 

e Ohmic, Schottky, and insulating contacts. 

e SRH, radiative, Auger, and surface recombination. 
e Impact ionization (local and non-local). 

e Floating gates. 

¢ Band-to-band and Fowler-Nordheim tunneling. 
¢ Hot carrier injection. 

¢ Quantum transport models 

e Thermionic emission currents. 


Table 6. Silvaco Atlas physical models (From [23, p. 


Cc. INPUT FILE STRUCTURE 


Silvaco Atlas receives input files through 





DeckBuild. 





The code entered in the input file calls Atlas to run with 





the following command: 


go atlas 


Following that command, the input file needs to follow 


a pattern. The command groups are listed in Figure 22. 
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Group Statements 


MESH 
REGION 
ELECTRODE 
DOPING 


1. Structure Specification 


MATERIAL 
MODELS 
CONTACT 
INTERFACE 


2. Material Models Specification 


3. Numerical Method Selection METHOD 


LOG 
SOLVE 
LOAD 
SAVE 


4. Solution Specification 


5. Results Analysis EXTRACT 
TONYPLOT 





Figure 22. Atlas command groups and primary statements (From 
[23% De 2-8 ])). 


Atlas follows the following format for statements and 


parameters: 





<STATEMENT> <PARAMETER>=<VALUE> 


The following line of code serves as an example. 




















DOPING UNIFORM N.TYPE CONCENTRATION=1.0e16 REGION=] ‘ 





OUTFILE=my.dop 








The statement is DOPING. The parameters are UNIFORM, 








N.TYPE, CONCENTRATION, REGION, and OUTFILE. There are four 











different type of parameters: real, integer, character, and 


logical. The back slash (\) serves the purpose of continuing 





the code in the next line. Parameters, such as UNIFORM, are 


30 


logical. Unless a TRUE or FALSE value is assigned, the 


parameter is assigned the default value. This value can be 





either TRUE or FALSE. The Silvaco Atlas manual needs to be 
referenced to identify the default value assigned to 


specific parameters. 


1. Structure Specification 





The structure specification is done by defining the 


mesh, the region, the electrodes and the doping levels. 


a. Mesh 


The mesh used for this thesis is two-dimensional. 
Therefore, only x and y parameters are defined. The mesh is 
a series of horizontal and vertical lines and spacing 
between them. From Figure 23, the mesh statements are 


specified. 


—- 


"i 
i 


x.mesh loc=-250 

x.mesh loc=-25 spac=2.5 
x.mesh loc=25 spac=2.5 
x.mesh loc=250 spac=25 


y-mesh loc=-2.9 spac=0.01 
y.mesh loc=-2.8 spac=0.03 
loc=-2.5 spac=0.05 
spac=0.05 
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ara 
-200 -100 0 100 200 
Microns 
Figure 23. Atlas mesh (From 2, p.18]). 
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The general format to define the mesh is: 


X.MESH LOCATION=<VALUE> SPACING=<VALUE> 








Y.MESH LOCATION=<VALUE> SPACING=<VALUE> 





For example, the x.mesh starting at -250 microns 


has spacing of 25 microns. That means it is relatively 





coarse. The x.mesh becomes finer between -25 and 25 microns 





with a spacing of 2.5 microns. The y.mesh is similarly 





defined. For example, at y.mesh of -2.9 microns, the spacing 
is 0.01 microns. Then at location y.mesh of -2.8 microns, 


the spacing changes to 0.03 microns. The mesh is coarser at 





y.-mesh location of —-1, when the spacing is 0.1. 


A coarse or fine mesh determines the accuracy of 





the simulation. A coarse mesh produces a faster simulation, 
but less accurate results. A fine mesh produces a slower 
simulation, but more accurate results. The areas that have a 
finer mesh, therefore, are of greatest interest in the 


simulation. 


b. Region 


After defining the mesh, it is necessary to define 





the regions. The format to define the regions is as follows: 








REGION number=<integer> <material_type> / 
<position parameters> 


From Figure 24, the code that defines the regions 





is identified. There are six regions defined. The limits of 








each region are explicitly identified in the x- and y-axis. 


The regions must then be given a material. 
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region n i : Bs v -8 y.-max=-2.5 


region n -mi X.max=-25 y.min=-2.9 y.max=-2.5 
region r x.min=25 x.max=250 y.min=-2.9 y.max=-2.5 
region n x.min=-250 x.max=250 y.min=-2.5 y.max=-2 
region x.min=-250 x.max=250 y.min=-2 y.-max=0 
region i a = y-min=-2.9 y.max=-2.8 














Figure 24, Atlas region (From [2, p. 19]). 


From Figure 25, the code defines the material for 





each region. Note that the color coding identifies the 





material. The regions have vertical and horizontal lines to 


mark their boundaries. 


region mate rial=InGaNas 


region mate rial=VYacuun 
region material =Vacuun 
region material =InGaNas 
region material =InGaNas 
region mate rial=Vacuun 





Materials 
InGaNAs 
Vacuum 
-200 -100 0 100 200 
Microns 
Figure 25. Atlas regions with materials defined (From [2, p. 
19]). 
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c. Electrodes 


The next structure specification corresponds to 


electrodes. Typically, in this simulation the only 








electrodes defined are the anode and the cathode. However, 
Silvaco Atlas has a limit of 50 electrodes that can be 


defined. The format to define electrodes is as follows: 








ELECTRODE NAME=<electrode name> <position_parameters> 


From Figure 26, the electrode statements are 
defined for the anode and the cathode. Note that the cathode 
is defined with gold as the material. The x and y dimensions 


correspond to region 6 previously defined. Meanwhile, the 





anode is defined at the bottom of the cell for the entire x 


range at y=0. 


electrode name=cathode material=Gold x.min=-25 x.max=25 Yy.min=-2.9 y.max=-2.8 


electrode name=anode x.mMin=-250 x.max=250 ¥.min=0 y.max=0 





-200 -100 Q 100 200 
Microns 


Figure 26. Atlas electrodes (From [2, p. 20]). 
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needs to be defined is doping. The format of the Atlas 


d. Doping 


The last aspect of structure specification that 


statement is as follows: 


DOP 








ING <distribution type> <dopant_type> / 


<position parameters> 


levels are defined. Doping can be n-type or p-type. The 


From Figure 27, the doping types and the doping 





distribution type can be uniform or Gaussian. 


Microns 


specification is next. From Figure 28, the materials model 


' 
~ 


‘ 
nv 











Net Doping (/cm3) 





doping uniform region=l n.type conc=2el8 
doping uniform region=4 n. type conc=1le20 
doping uniform region=5 p.type conc=d8elé8 


~200 ~100 0 100 200 
Microns 


Figure 27. Atlas doping (From [2, p. 21]). 


2. Materials Model Specification 


After the structure specification, the materials model 











specification is broken down into material, models, contact, 





and interface. 
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MATERIAL 


: Se MODELS 
2. Material Models Specification ——_—_—_ CONTACT 
INTERFACE 
Figure 28. Atlas material models specification (After [23, p. 


2-8]). 


a. Material 


The format for the material statement is as 


follows: 


MATERIAL <localization> <material_definition> 











Below are three examples of the material 


statement: 


MATERIAL MATERIAL=Silicon EG300=1.1 MUN=1200 














MATERIAL REGION=4 TAUNO=3e-7 TAUPO=2e-5 








MATERIAL NAME=base NC300=4e18 








In all examples, when MATERIAL appears first, it 





is considered the statement. When MATERIAL appears a second 


time in the first example, it is considered a localization 





parameter. In the second and third examples, the 





localization parameters are REGION and NAME, respectively. 
Various other parameters can be defined with the material 
statement. Examples of these parameters are the band gap at 
room temperature (EG300), electron mobility (MUN), electron 


(TAUNO) and hole (TAUPO) recombination lifetimes, conduction 





band density at room temperature (NC300), among others. 


b. Models 





The physical models fall into five categories: 
mobility, recombination, carrier statistics, impact 
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ionization, and tunneling. The syntax of the model statement 


is as follows: 





MODELS <model flag> <general parameter> / 
<model dependent parameters> 


The choice of model depends on the materials 


chosen for simulation. 


The example below activates several models. 





MODELS CONMOB FLDMOB SRH 








CONMOB is the concentration dependent model. 
FLDMOB is the parallel electric field dependence model. SRH 














is the Shockley-Read-Hall model. 


c. Contact 


Contact determines the attributes of the 


electrode. The syntax for contact is as follows: 





CONTACT NUMBER=<n> |NAMF=<ename>|ALL 


The following is an example of the contact 


statement. 


CONTACT NAMF=anode current 


d. Interface 


The semiconductor or insulator boundaries are 








determined with the interface statement. The syntax is as 


follows: 


INTERFACE [<parameters>] 





The following example shows the usage of the 





interface statement. 
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NTERFACE X.MIN=—-4 X.MAX=4 Y.MIN=-0.5 Y.MAX=4 \ 








OF=1lel0 S.N=le4 S.P=1le4 


The max and min values determine the boundaries. 








The QF value specifies the fixed oxide charge density (cm 
ay The S.N value specifies the electron surface 
recombination velocity. S.P is similar to S.N, but for 


holes. 
3. Numerical Method Selection 


After the materials model specification, the numerical 





method selection must be specified. From Figure 29, the only 
statement that applies to numerical method selection is 


method. 


3. Numerical Method Selection —_—_— METHOD 


Figure 29. Atlas numerical method selection (After [23, p. 2- 
8.])) 


There are various numerical methods to calculate 





solutions to semiconductor device problems. There are three 


types of solution techniques used in Silvaco Atlas: 
e decoupled (GUMMEL) 
e fully coupled (NEWTON) 


e BLOCK 


The GUMMEL method solves for each unknowns by keeping 


all other unknowns constant. The process is repeated until 





there is a stable solution. The NEWTON method solves all 
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unknowns simultaneously. The BLOCK method solves some 
equations with the GUMMEL method and some with the NEWTON 
method. 


The GUMMEL method is used for a system of equations 





that are weakly coupled and there is linear convergence. The 
NEWTON method is used when equations are strongly coupled 


and there is quadratic convergence. 


The following example shows the use of the method 


statement. 


METHOD GUMMEL NEWTON 








In this example, the equations are solved with the 





GUMMEL method. If convergence is not achieved, then the 


equations are solved using the NEWTON method. 


4. Solution Specification 


After completing the numerical method selection, the 








solution specification is next. Solution specification is 
broken down into log, solve, load, and save statements, as 


shown in Figure 30. 


LOG 
4. Solution Specification — rig 
SAVE 
Figure 30. Atlas solution specification (After [23, p. 2-8]). 
a. Log 


LOG saves all terminal characteristics to a file. 





DC, transient, or AC data generated by a SOLVE statement 


after a LOG statement is saved. 
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The following shows an example of the LOG 


statement. 


LOG OUTFILE=myoutputfile.log 








The example saves the current-voltage information 


into myoutputfile.log. 


b. Solve 


The SOLVE statement follows the LOG statement. 
SOLVE performs a solution for one or more bias points. The 


following is an example of the SOLVE statement. 


SOLVE B1=10 B3=5 BEAM=1 SS.PHOT SS.LIGHT=0.01 \ 








MULT.F FREQUENCY=le3 FSTEP=10 NFSTEP=6 


Bl and B3 specify the optical spot power 


associated with the optical beam numbers 1 and 3, 





respectively. The beam number is an integer between 1 and 


10. BEAM is the beam number of the optical beam during AC 








photogeneration analysis. SS.PHOT is the small signal AC 








analysis. SS.LIGHT is the intensity of the small signal part 
of the optical beam during signal AC photogeneration 


analysis. MULT.F is the frequency to be multiplied by FSTEP. 





NFSTEPS is the number of times that the frequency is 


incremented by FSTEP. 


c. Load and Save 








The LOAD statement enters previous solutions from 


files as initial guess to other bias points. The SAVE 








statement enters all node point information into an output 





file. 
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The following are examples of LOAD and _ SAVE 


statements. 


SAVE OUTF=SOL.STR 





In this case, SOL.STR has information saved after 
a SOLVE statement. Then, in a different simulation, SOL.STR 


can be loaded as follows: 





LOAD INFILE=SOL.STR 











5. Results Analysis 


Once a solution has been found for a semiconductor 
device problem, the information can be displayed graphically 
with TonyPlot. Additionally, device parameters can be 


extracted with the EXTRACT statement, as shown in Figure 31. 


5. Results Analysis we 2 EX RACT 
TONYPLOT 
Figure 31. Atlas results analysis (After [23, p. 2-8]). 


In the example below, the EXTRACT statement obtains the 





current and voltage characteristics of a solar cell. This 





information is saved into the IVcurve.dat file. Then, 





TonyPlot plots the information in the IVcurve.dat file. 








EXTRACT NAME="iv" curve(v."anode", I."cathode") / 
OUTFILE="IVcurve.dat" 














TONYPLOT IVcurve.dat 





Figure 32 shows the sample IV curve plotted by 
TonyPlot. 
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=) Tonyplot 
File Edit Plot Tools Production Help 


SMWASBS LOAeA AQ es 





ATLAS 
Data from SingleCell_webf log 


*——* Cathode Current (4) 


oot 006 ons O41 
Anode Voltage (¥) 


Tonyplot 3.4.4.R © Silvaco 2007 








Figure 32. Sample TonyPlot IV curve. 


D. CONCLUSION 


This chapter presented an introduction to Silvaco 
Atlas, the structure of the input files, and some of its 


statements. With these basic tools, the next chapter 





introduces some basic information on wurtzite Indium Gallium 


Nitride before proceeding to the thesis simulation. 
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IV. INDIUM GALLIUM NITRIDE 


A. A FULL SPECTRUM PHOTOVOLTAIC MATERIAL 





Prior to 2001, it was thought that wurtzite Indium 
Nitride had a band gap of approximately 1.9 eV. In 2001, it 








was discovered that Indium Nitride had a much smaller band 











gap. Davydov, et al., concluded in [24] that Indium Nitride 





had a band gap of approximately 0.9 eV. Additionally, 








Davydov, et al., presented in [25] that wurtzite Indium 


Gallium Nitride band gaps in the range x=0.36 to x=l1 





supported findings in [24]. The band gaps for Indium Gallium 











Nitride, with Indium Nitride concentrations ranging from 0.0 


to 1.0 are presented in Figure 33. 


Energy (eV) 


- M.H.Kim 

- K.P.O'Donnell 

- A.Klochikhin 

- V.Yu.Davydov (2002) 











0 
0.0 0.2 0.4 0.6 0.8 1.0 
INN fraction, x 
Figure 33. InGaN band gap as a function of In composition 
(After [25]). 
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Later on, research conducted at Lawrence Berkeley 





National Laboratory presented that the Indium Nitride 
fundamental band gap was approximately 0.77 eV at room 
temperature. Since Gallium Nitride has a band gap of 


approximately 3.4 eV at room temperature, then Indium 





Gallium Nitride can have a band gap ranging from 0.77 eV to 








3.4 eV by changing the percent composition of Indium and 








Gallium within Indium Gallium Nitride. Figure 34 confirms 


that Indium Gallium Nitride follows the pattern of ranging 








from 0.77 eV to 3.4 eV. 



















4} +4 
| @ LANL, 2002 
3 | @ Literature 13 
> a 
2 e g 
Co] MO. ee | 2% 
en = 
ArT Ss 
1 |} 10 
5 4 3 2 1 0.0 0.2 0.4 0.6 0.8 1.0 
Solar Flux Ga content in 
(1071 photons/secim/m) In,, Ga,N alloy 
Figure 34. InGaN band gap and solar spectrum comparison 





(After [26]). 


Figure 34 also presents the solar spectrum on the left- 





hand side. Therefore, it is clear that the visible spectrum 





has a near-perfect match with the Indium Gallium Nitride 





range of band gaps. 





It should be noted that research conducted in [27] 





presents a Indium Nitride band gap of 1.740.2 eV. It is 








stated in [27], “These investigations clarify that a band 


44 





transition around 1.7 eV does exist in wurtzite InN although 


in this energy range optical material responses are hardly 


ever found. 
fundamental 


concluded that 














If the 1.7 eV band transition is indeed the 








bandgap in InN at room temperature it is 





defect bands, grain boundaries, dislocations 





and/or the very conductive surface are contributing to the 








lower energy optical responses in InN. The identification of 


the various 


investigations.” 





In 


follow 


up 





defects in InN will be subject of future 


research, [28] explains the possible 





reason for the higher Indium Nitride band gap: “Owing to 








InN’s exceptional propensity for n-type activity, in many 


instances the larger apparent values of the band gap can be 





attributed to 





the Burstein-Moss shift of the optical 


absorption edge resulting from the occupation of conduction 





band states by free electrons.” 


Figure 35 shows evidence that the fundamental band gap 


of indium n 


itride 





is approximately 0.7 eV, and not near 2 





eV. The opt 








tical absorption (blue curve) has an onset at 

















approximately 0.8 eV. There is no change in absorption 
around the 2 eV range. The photo luminescenc (red curve) 
shows a peak at the band edge of approximately 0.8 eV. 
Finally, the photo-modulated reflectance (green curve) 
presents a transition close to 0.8 eV. These three curves 
support a band gap for indium nitride of approximately 0.8 
eV. The 1.9-2.0 eV band gap is not evident in Figure 35. 
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Figure 35. Evidence of 0.7 eV band gap for indium nitride 
(From [29]). 


The controversy about indium nitride’s fundamental band 





gap is not the subject of this thesis. This research 
proceeded accepting the fundamental band gap of 0.77 eV for 


indium nitride. 





According to [30], the following formula provides an 








approximation of Indium Gallium Nitride band gap: 


Ec(x) =3.42x+0.77(1— x) —1.43x(1— x) 





Where Ee is the InGaN band gap, 3.42 eV is the GaN band 





gap, 0.77 eV is the InN band gap, 1.43 eV is the bowing 


parameter b, x is the Ga concentration, and (1-x) is the In 





concentration. 


When plotting the band gap formula for InGaN using a 





Matlab script, the Figure 36 is obtained. 
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Indium concentration vs InGaN Band Gap 
3.5 


N 
in 


InGaN Band Gap (eV) 
in N 


0.5 
oO o1 O2 O38 O4 O58 O68 OF o8 og 1 
Indium concentration (unitless ratio) 
Figure 36. InGaN band gap as a function of In concentration. 








The band gaps used in the Silvaco Atlas simulations 


were found by modifying the Matlab script for Figure 35. 


B. RADIATION-HARD SEMICONDUCTOR MATERIAL 


Radiation-hard materials are essential for space 


applications. When compared to other Gallium Arsenide and 








Indium Gallium Phosphide, Indium Gallium Nitride is able to 





withstand a greater amount of radiation. According to [31], 





InGaN “retains its optoelectronic properties at radiation 


damage doses at least two orders of magnitude higher than 





the damage thresholds of the materials (GaAs and GalInP) 
currently used in high efficiency MJ cells”. Radiation 


levels in [31] were 1 MeV electron, 2 MeV proton, and 2 Mev 








alpha particle irradiation. 
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This indicates that InGaN is not only potentially a 


high-efficiency photovoltaic material, but it is also able 





to withstand the harsh space environment. 


Cc. INDIUM GALLIUM NITRIDE CHALLENGES 


There are some significant issues that currently 











prevent InGaN from being used as photovoltaic material. As 
of this writing, no PN junction has been built. The PN 
junction is essential to create photodiodes. Indium Nitride 








has the highest electron affinity among all semiconductors 








[32]. Hence, Indium Nitride has a tendency to be n-type. It 








is much more difficult to build p-type material with Indium 














Nitride. However, videnc of p-type doping in Indium 











Nitride and Indium Gallium Nitride has been reported in [33] 





using Magnesium. This is a step in the right direction. 


Further progress in this area dictates the creation of 





Indium Gallium Nitride PN junctions. 











Having reviewed the characteristics of Indium Gallium 





Nitride, the next chapter covers the simulation of Indium 


Gallium Nitride using Silvaco Atlas. 
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V. SIMULATION OF INDIUM GALLIUM NITRIDE IN SILVACO 
ATLAS 


One of the primary motivations behind this research was 





to find a way to make solar cells significantly more 
efficient. When researching high-efficiency solar cells, the 


improvements were typically incremental over previous 





results. Solar cells were first designed as single-junction. 


The most popular photovoltaic material has been Silicon due 





to its abundance and relatively low cost. The maximum 








efficiency of polycrystalline Silicon is approximately 





19.8%. Commercial silicon solar cells have a maximum 
efficiency of approximately 13%. Improvements in efficiency 








can be obtained in single-junction solar cells by using 


Gallium Arsenide. Monocrystalline Gallium Arsenide has a 





maximum efficiency of approximately 25.1% [10, p. 195]. 





By designing dual-junction solar cells, the efficiency 
can again be improved. A dual-junction cell with Gallium 
Phosphide as the top junction and Gallium Arsenide as the 


bottom cell can have an efficiency of 30.3% [10, p 302]. 





Increasing the number of junctions makes the cells more 
complex in design and the gains begin to diminish. Triple- 


junction cells have efficiencies of approximately 31%. Quad- 





























junction cells have efficiencies of approximately 33%. It is 
expected that five- and six-junction solar cells can have 


efficiencies of approximately 35% [36]. 





This simulation uses Indium Gallium Nitride as 


photovoltaic material in a simulation. With its full- 





spectrum characteristics, the efficiency of Indium Gallium 


Nitride has the potential to surpass previous results. 
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A. SINGLE-JUNCTION SOLAR CELL 


Figure 37 shows a graphical representation of the 
Single junction solar cell. The thickness of the emitter is 


0.01 wm and the doping is 1x10'°/cm®. The thickness of the 





base is 3 wm and the doping is 1x10*°/cm’. 


Emitter ntInGaN O.0lyzm  lel6cm™ 


Base ptInGaN 3m 








Figure 37. Simple single-junction InGaN solar cell. 





The first step is to develop single junction solar 


cells. Optical data in the form of dielectric constants was 





found in [34] for Ino.20Gao.g0N. This corresponds to a 


calculated band gap of 2.66 eV. Additional optical data was 








found in [35] for IMno.57Gao.43N (calculated Eg=1.60 eV), 








Ino. 6gGao.32N (calculated Eg=1.31 ev), and Ino.7gGao.22N 
(calculated Eg=1.1 eV). The dielectric constants were then 
converted to index of refraction (n) and extinction 


coefficient (k) using a Matlab script (Appendix B). The band 
gap and the optical data were entered into the input deck. 
The simulation ran in Silvaco Atlas. From the log file, the 


current and voltage data was extracted. 





In the single-junction case, each of the cells had AMO 





as the input spectrum. Figure 38 shows the IV curves of four 





single-junction solar cells. The IV curve pattern 


corresponds to empirical findings of other multijunction 
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solar cells. Typically, the cell with the highest voltage 
has the lowest current. The cell with the lowest voltage has 


the highest current. 


Solar Cell I'Y¥ Curve 


0.025 





0.02 
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= 
o 
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= 
7) 
= 
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a 
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0 0.5 1 1.5 2 25 
Voltage (V) 
Figure 38. Four single-junction IV curves. 








Table 7 shows Isc, Voc, fill factor, and efficiency of 





the single-junction solar cells presented in Fig. 36. 



























































Material Isc Voc (V) Fill Efficiency 
(mA/cm’) factor 
Ino.20Gao.s80N 17.613 Diet 83.72% 24.32% 
Ino.57Gao.43N | 18.395 Le eS 78.80% 12.35% 
Ino.6gGao0.32N 20.453 0.91 76.38% 10.52% 
Ino.7gGao.22N 20.848 0.72 73.73% 08.19% 
Table 7. Efficiencies of four single-junction InGaN cells. 
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The efficiencies shown in Table 7 indicate that a high- 





efficiency, single-junction cell with Ino.290Gao.gpN has good 
prospects. A 24.32% efficiency compares favorably to the 


single-junction Gallium Arsenide efficiency of 25.1%. The 





InGaN cells simulated can be optimized by changing the 
thickness of the junction and changing the doping levels. 
Adding a window, a back surface field (BSF) and a buffer can 
also help in improving the efficiency results. These 


improvements can be the subject of future research. 


B. DUAL-JUNCTION SOLAR CELL 


Figure 39 shows a graphical representation of a dual- 








junction InGaN solar cell. The thicknesses and doping levels 


remain identical to the single-junction case. 





Emitter notInGaN 0.01pm  lelé6cm™ 


Base ptInGaN 3m lel6cm™ Eg=2.66 eV 


Emitter ntInGaN 0.01pm  lel6cm™ 


Base ptInGaN 3m lel6cm™ 








Figure 39. Simple dual-junction InGaN solar cell. 


From the discussion in the single-junction solar cell 


section, it should be clear that the best combination of 





InGaN band gaps is 2.66 eV and 1.60 eV, since they have the 
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highest efficiencies. It should be noted that for this 


Simulation, the AMO spectrum was provided for the top 





junction Ino.20Gao.goN (Eg=2.66 eV). The bottom junction 











Ino.57Gag.43N (Eg=1.60 eV) received the AMO spectrum minus the 


spectrum absorbed by the top cell. 





Figure 40 shows the IV curves of the dual-junction 


solar cell. Note that compared to Figure 38, the IV curve of 











Eg=2.66 eV stayed the same, while the IV curve of Eg=1.60 eV 
had a current drop. This current decrease was expected 


because the input spectrum was diminished by the top cell. 


Solar Cell I'Y¥ Curve 
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Figure 40. Dual-junction InGaN solar cell IV curve. 











Since the junctions are in series, the overall IV curve 








is limited in its current level by the IV curve with the 





lowest current. The voltages are added accordingly. 


Table 8 shows the efficiency of the dual-junction solar 


cell. 
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Material Isc Voc (V) Fill Efficiency 
(mA/cm?) factor 
Ino.20Gao.goN Egy ae seed 3.38 83.25% 35.58% 
Ino,57Gao,43N 
Dual 
junction 
Table 8. Dual-junction InGaN efficiency. 

There is a Significant increase in efficiency from 

single-junction (24.32%) to dual junction (35.58%) InGaN 





solar cells. This improvement alone should be enough to 








continue to pursue the development of InGaN solar cells. 


Searching for further improvements, three- and four- 


junction solar cells are examined next. 





Cc. THREE-JUNCTION SOLAR CELL 


Figure 41 shows a graphical representation of a three- 





junction InGaN solar cell. Thicknesses and doping levels 





stayed the same as the single- and dual-junction cases. From 


the single-junction solar cell discussion, the three most 





efficient cells wer selected (Eg=2.66 eV, Eg=1.60 eV, 





Eg=1.31 eV). The top cell Inpo.20Gao.goN (Eg=2.66 eV) received 








the AMO spectrum. The cell Ino.57Gao.43N (Eg=1.60 eV) received 








the AMO spectrum minus the spectrum absorbed by the top 











cell. The cell Ino.6sGao.32N (Eg=1.31 V) received the AMO 


spectrum minus the spectrum absorbed by the top two cells. 








Figure 42 shows the IV curve of the three-junction 


solar cell. Note that the current level of the top two cells 








is identical to the dual-junction case. However, the current 
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level of the bottom cell is lower than the single junction 


case for the same concentration. Similarly to the dual- 





junction case, the current of the overall IV curve is 


limited by the individual junction with the lowest IV curve. 





The voltages are added accordingly, since the three 


junctions are in series. 


Emitter ntInGaN OOlgm  lel6cm™ 


Base  ptInGaN 34m lel6cm™ Eg=2.66 eV 


Emitter ntInGaN OOlgm  lel6cm™ 


Base ptInGaN 3ym lel6cm™ Eg=1.60 eV 


Emitter ntInGaN O0Olgm  lel6cm™ 


Base ptInGaN 3m lel6cm™ Eg=1.31 eV 








Figure 41. Simple three-junction InGaN solar cell. 


55 





Current Density (A¥cm?) 























Solar Cell I'¥ Curve 


—— Eg=2.66 eV 
——— Eg=1.60 eV 
—— Eg=1.31 eV 
Overall 
© Max power 
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Voltage (V) 
Figure 42. Three-junction InGaN solar cell IV curve. 
Table 9 shows the efficiency of the three-junction 
InGaN solar cell. 
Material Isc Voc (V) Fill Efficiency 
(mA/ cm’) factor 
Ino.20Gao.g0N 14.12 4.27 87.23% 38.90% 
Ino.57Gao.43N 
Ino. 6sGao.32N 
Three 
junction 
Table 9. Three-junction InGaN efficiency 








The increase in efficiency from dual-junction to three- 


junction is 


35.58% 


to 


38.90%. The 


efficiency can be attributed to the drop in 
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modest 





increase in 


Isc from 17.11 





mA/cm* to 14.12 mA/cm’. However, the increase in Voc from 





3.38 V to 4.27 ensured that the efficiency would increase. 





The last case to be examined is the quad-junction InGaN 





solar cell. 


D. QUAD-JUNCTION SOLAR CELL 





Figure 43 shows a graphical representation of a quad- 





junction InGaN solar cell. Thicknesses and doping levels 


remained the same as the single-, dual-, and three-junction 





cases. Spectrum input followed the pattern of the dual- and 
three-junction cases. As expected, the current dropped for 


all cases where AMO had been reduced by the top cells. 








Figure 44 shows the IV curve for a quad-junction InGaN 





solar cell. Note that the IV-curves for the top three cells 











are identical to the three-junction case. The bottom cell 





has a lower current level compared to the single junction 


Case. 


Table 10 shows the efficiency of the gquad-junction 








InGaN solar cell. 


Dh 


ptInGaN 34m lel6cm™ Eg=2.66 eV 


ptInGaN 3um lel6cm™ Eg=1.60 eV 


ptInGaN 34m lel6cm™> Eg=1.31 eV 


ptInGaN 34m lel6cm™ Eg=1.10 eV 





Figure 43. Simple quad-junction InGaN solar cell 
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——— Eg=1.10 eV 
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Voltage (Vv) 
Figure 44. Quad-junction InGaN solar cell IV curve 
Material Isc Voc (V) Fill Efficiency 
(mA/cm?) factor 
Ino.20Gao.s0N 12.88 4.98 87.86% 41.69% 








5 


0.57GA0.43N 


Ino. 6sGao.32N 














Ino.7gGao.22N 


Quad-— 


junction 

















Table 10. 


For completeness, 


below for the quad-junction case. 


Quad junction 





performed with the help of a Matlab script. 
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InGaN efficiency 


efficiency calculations are provided 


These calculations were 





A (4.4736 V)=0.0564 








Pmax=Imp Vmp=(0 0126 

















ine cm? 

. IV 0.0564 
FF= max. _.Amp-¥ mp a cm. =87.86% 

Isc Voc Isc Voc (0.0129 5 )(4.9793V) 

cm 

Poe Teevy, 0.0564 > 
n= ae — +mp , a = WW =0.4169=41.69% 

= Pe Gage 

cm 


In order to get a perspective on how this result 





compares to actual solar cells in production, Figure 45 is 


presented. 

















35% 
=@— Si Space PV 
—H- GaAs SJ nJ 
xTJ 45 to 6J 
it~ GalnP/GaAs/Ge MJ GainP/GaAsiGe 
30% A 
A XTJ 
> Ow TJ uTJ 
= a ag 
2 25% 
FS am DJ TJ 
Ww o GalnP/GaAs/Ge GalnP/GaAs/Ge 
£2 sJ 
3 z 20% GaAsiGe 
SJ 
oS GaAsiGaAs 
— — 
3 = BSF - BSR Si 
CE 15% 
o ad Textured Si Thin Si (84 um) 
20 
go J 
SS Large Area Si 
£g 10% {2 x 6 cm) fs 
ry GalnP/GaAs/Ge MJ 
a : 
ao — Generation 
: . rati 
=o GaAs SJ Generation 
0% 
1960 1965 1970 1975 1980 1985 1990 1995 2000 # 2005 2010 2015 
Date of First Flight 
Figure 45. Spectrolab’s solar cell efficiencies (From [34]). 
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The label nJ stands for new generation four to six- 


junction cells. XTJ is a designation for a triple junction 





cell. UTJ stands for Ultra Triple Junction. ITJ means 











Improved Triple Junction. TJ is Triple Junction. DJ is Dual 





Junction. SJ is single junction. 


Figure 45 shows the efficiency progression of a leading 
solar cell manufacturer (Spectrolab). The first step 
consists of single junction silicon solar cells. These cells 
were in production from the mid-1960s until the early 1990s. 
Their efficiency was close to 15%. The next step consisted 


of single junction gallium arsenide solar cells. These cells 











had efficiencies between 15% and 20%. In the mid-1990s, dual 








junction solar cells appeared. Efficiencies wer slightly 
above 20%. Triple junction cells appeared in the early 2000s 
with efficiencies nearing 30%. The projections are to 
produce quad-junction solar cells within the next decade 


with efficiencies just under 35%. 


Further comparisons can be made with another 


simulation. Modeling in [5] used physics equations to 





calculate Isc and Voc. The input spectrum used in that 





simulation was AM1.5 instead of AMO. The InGaN band gap 
formula used was from [37]. This thesis used the band gap 
formula from [30]. Figure 46 shows a comparison of the two 


formulas. Both formulas seek to duplicate findings from 








actual band gap measurements. Therefore, both formulas are 


approximations only. 


The simulation from [5] obtained the results presented 





in Table 11. The highest efficiency obtained in that 





simulation is 40.346% with a six-junction InGaN solar cell. 
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Indium concentration vs InGaN Band Gap 


Ref 25 forrnula (used in this thesis) 
Ref. 33 formula (used in Ref. 32 simulation) 


InGaN Band Gap (eV) 





“0 o1 a2 03 O4 05 O68 OF O8 O9 1 
Indium concentration (unitless ratio) 








Figure 46. Comparison of InGaN band gap formulas 
Number of Toe (A/cm?) Voc (V) Output peak power Efficiency (%) 
junctions in the cell (mW/cm?) 

2 27 1.22 26.48 27.485 
3 18.2 2.22 32.42 33.642 
4 13.5 3.28 35.48 36.825 
5 11 4.28 37.68 39.105 
6 9.1 5.34 38.88 40.346 

Table 11. InGaN efficiency results (From [5]). 





The differences in efficiency results can be partly 


attributed to the difference in band gap formula used. For 





example, IMno.20Gao.goN has a band gap of 2.66 eV with the 








formula used in this thesis. The formula from [37] yields a 





band gap of 1.897 eV. If the Ino.20Gao.goN band gap of 1.897 





eV is substituted in the single-junction solar cell 


simulation of this thesis, the efficiency decreases 








significantly. Figure 47 compares the IV output for 





Ino.20Gao.goN at Eg=1.897 eV and Eg=2.66 ev. 
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Solar Cell I'Y¥ Curve 


0.025 











0.02 
xo 
= 
oOo 
z 0.015 
2 
a 
= 
ao 
e 
€ O01 
® 
5 
oOo 
0.005 
0 
0 0.5 1 1.5 2 2.5 
Voltage (V) 
Figure 47. IV curve for Ino.20Gao.goN using different band gaps 





The efficiency of the Ino.20Gao.goN (2.66 eV) single- 
junction solar cell simulated in this thesis was 24.32%. 
Changing the band gap to 1.897 eV yields an efficiency of 
14.94%. This is a significant decrease. However, the band 
gap for the other three junctions increases with the formula 


used in [37]. The quad-junction simulation was run with the 








new band gaps. Figure 48 shows the IV curve for the quad- 
junction solar cell. The efficiency for the quad-junction 


solar cell increases to 43.62%. It is difficult to state 








which band gap model is more correct. 
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Solar Cell I'Y¥ Curve 


—— Eg=1.897 
——— Eg=1.825 
—— Eg=1.69 
—— Eg=1.43 
Overall 
© Max paint 


Current Density (A¥cm?) 











0 05 1 #15 2 25 3 35 4 45 #5 565 
Voltage (Vv) 
Figure 48. Quad-junction InGaN solar cell IV curve using 


calculated band gaps from [37] formula. 











Material Isc Voc (V) Fill Efficiency 
(mA/cm?) factor 
Ino.20Gao.goN | 12.9 5.194 88.06% 43.62 


Ino.57Gao.43N 


Ino.6sGao.32N 














Ino.7gGao.22N 


Quad-— 


junction 




















Table 12. Quad junction InGaN efficiency using band gaps 
from [37] calculations. 
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A common finding in [5] and this thesis is that both 


simulations show that 





Indium Gallium Nitride multijunction 





solar cells can provide a significant improvement in solar 





cell efficiency. However, this 





Technology Computer Aided Desi 


Atlas, has been used 


The 





Indium 


results of this thesis simulation demonstrate that 


to simulate 


is the first time that a 





gn (TCAD), such as Silvaco 





InGaN solar cells. 








Gallium Nitride is 


potentially an excellent 


semiconductor photovoltaic material. 





The material science 


research in the future can confirm these outcomes in the 


future. 


Product 


the firs 


are demonstrated with actual 


complexi 


Tunnel j 








ty of dual-junction sol 








to be explored as well. 
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Indium 





unctions made up of Indium Gall 


lar cel 


tion of single-junction solar cells appears to be 





t step to be taken. Once the photovoltaic properties 


Gallium Nitride, the 


ls can be addressed. 











ium Nitride may need 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. RESULTS AND CONCLUSIONS 











Indium Gallium Nitride is a semiconductor material with 
potential to be used in photovoltaic devices. A new 
simulation was performed using Silvaco Atlas. The results of 








the quad-junction Indium Gallium Nitride solar cell indicate 





that a new high-efficiency material should be produced. 


Current state-of-the-art multijunction solar cells have 








efficiencies in the 30-33% range. The 41% efficiency 
predicted in this thesis is the highest of all simulations 





performed at the Naval Postgraduate School using Silvaco 





Atlas. 


B. RECOMMENDATIONS FOR FUTURE RESEARCH 


There are multiple areas that can be explored in future 





research. This thesis focused on the development of a solar 








cell model that emphasizes the use of optical constants 








(refraction and extinction coefficients n and k). The model 
used default settings for other parameters, such as 
permittivity, affinity, radiative recombination rate, 


electron and hol lifetimes, lectron and hole density of 








states, and lattice constants. One area of future research 


is to obtain measured data for the above parameters. This 





subject is critical in improving the model. The other option 





is to interpolate from known values. Additionally, the 


inclusion of tunnel junctions in the simulation is of great 





value. 
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Other areas of research include finding different 














semiconductor materials that offer potential for high- 
efficiency solar cells. Lawrenc Berkeley National 
Laboratory is also working with multiband materials, such as 


Zinc Manganese Tellerium Oxide (ZnMnTeO). This type of 
material offers a single junction solar cell, but there are 


multiple band gaps within the single junction. 








One more recommended path is to optimize the physical 





parameters of the solar cell. The thickness of each of the 


layers can be optimized to produce a better efficiency. The 





band gap distribution for a quad-junction solar cell can be 





improved to obtain higher efficiency. It should be noted 





that obtaining the correct optical data for these band gaps 


is essential. 
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APPENDIX A: SILVACO ATLAS INPUT DECK 





The code for the quad-junction InGaN solar cell was 


broken down into four single-junction solar cells. This was 








don becaus the material properties of InGaN had to be 


changed for each of the junctions. 


Original code for the Silvaco Atlas input deck was 
obtained from Michalopolous [1], Bates, [2], Green [3], and 


Canfield [4]. Modifications to the code were made to support 





this thesis. 


A. TOP JUNCTION: INo.20GAo.go0N, EG=2.66 EV 


go atlas 

set cellWidth=5.000000e+002 

set capWidthpercent=8.000000e+000 
set divs=1.000000e+001 

set contThick=1.000000e-001 


set capDop=1.000000e+020 
set windowThick=0.01 
set winDop=2.15e17 


set capThick=3.000000e-001 


mitterThick=0.01 

#changed emitDop from lel6 to 1e20 
set emitDop=1lel16 

set baseThick=3.19467 

#changed basDop from lel6é to 1e20 

















set baseDop=1lel16 

set bsfThick=0.03533 

set bsfDop=2.15e19 

set cellWidthDiv=$cellWidth/Sdivs 
set width3d=100e6/$cellWidth 

set capWidth=0.01*ScapWidthpercent*S$cellWidth/2 
set capWidthDiv=ScapWidth/ (Sdivs/2) 
set cellWidthHalf=$cellWidth/2 

set bsfLo=0 

set bsfHi=S$bsfLo-SbsfThick 

set bsfDiv=SbsfThick/Sdivs 

set baseLo=SbsfHi 

set baseHi=S$baseLo-SbaseThick 

set baseMid=SbaseLo-SbaseThick/2 
set baseDiv=SbaseThick/Sdivs 

set emitterLo=SbaseHi 
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set emitterHi=SemitterLo-SemitterThick 
set emitterDiv=SemitterThick/S$divs 

set windowLo=SemitterHi 

set windowHi=SwindowLo-SwindowThick 
set windowDiv=$windowThick/Sdivs 

set capLo=$windowHi 

set capHi=ScapLo-ScapThick 


#set capDiv=ScapThick/S$divs 


set contLo=$capHi 
set contHi=$contLo-ScontThick 
set contDiv=ScontThick/Sdivs 





set lightY=SemitterHi-5 


mesh width=S$width3d 

## X-Mesh 

x.mesh loc=-$cellWidthHalf spac=$cellWidthDiv 
-mesh loc=-ScapWidth spac=ScapWidthDiv 

-mesh loc=ScapWidth spac=S$capWidthDiv 

-mesh loc=$cellWidthHalf spac=ScellWidthDiv 


x x x 


# Y—-Mesh 

Top. contact 

y-mesh loc=ScontHi spac=0 
y-mesh loc=ScontLo spac=0 
Cap 
Window 

y-mesh loc=SwindowHi spac=SwindowDiv 
y-mesh loc=SwindowLo spac=SwindowDiv 
Emitter 

y.mesh loc=SemitterLo spac=SemitterDiv 
Base 
y-.mesh loc=SbaseMid spac=SbaseDiv 
BSF 
y.mesh loc=SbsfHi spac=SbsfDiv 
y.mesh loc=SbsfLo spac=SbsfDiv 

















## Regions 

# Cap 

#region num=8 material=Vacuum x.min=—-ScapWidth x.max=ScapWidth y.min=ScontHi 
y.max=ScontLo 

region num=1 material=InGaN x.min=-ScapWidth x.max=ScapWidth y.min=ScapHi 
y.max=ScapLo x.comp=0.20 

region num=2 material=Vacuum x.min=-S$cellWidthHalf x.max=—-ScapWidth 
y.min=ScontHi y.max=ScapLo 

region num=3 material=Vacuum x.min=ScapWidth x.max=ScellWidthHalf y.min=ScontHi 
y.max=ScapLo 

# Window [for Ge cell, use AlGaAs with x.comp=0.2] 

region num=4 material=InGaN x.min=-ScellWidthHalf x.max=S$cellWidthHalf 
y.-min=SwindowHi y.max=SwindowLo x.comp=0.20 


#region num=4 material=AlGaAs x.comp=0.2 x.min=—-S$cellWidthHalf 
x.max=ScellWidthHalf y.min=SwindowHi y.max=SwindowLo 

# Emitter 

region num=5 material=InGaN x.min=-$cellWidthHalf x.max=ScellWidthHalf 
y-min=SemitterHi y.max=SemitterLo x.comp=0.20 

# Base 
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region num=6 material=InGaN x.min=-$cellWidthHalf x.max=ScellWidthHalf 
y.min=SbaseHi y.max=SbaseLo x.comp=0.20 

# BSF 

region num=7 material=InGaN x.min=—-S$cellWidthHalf x.max=ScellWidthHalf 
y-.min=SbsfHi y.max=SbsfLo x.comp=0.20 


## Electrodes [for InGaP cell, add cathode (gold) and remove cathode 
(conductor) ] 





electrode name=cathode material=Gold x.min=—-ScapWidth x.max=ScapWidth 
y.min=ScontHi y.max=ScontLo 
#electrode name=cathode x.min=—-S$cellWidthHalf x.max=ScellWidthHalf 


y.min=SwindowHi y.max=SwindowHi 
lectrod name=anod x.min=-ScellWidthHalf x.max=S$cellWidthHalf y.min=SbsfLo 
y.max=SbsfLo 














# Doping [for InGaP cell, uncomment cap doping] 
Cap 

doping uniform region=1 n.type conc=ScapDop 
Window 

doping uniform region=4 n.type conc=$winDop 
Emitter 

doping uniform region=5 n.type conc=SemitDop 
Base 

doping uniform region=6 p.type conc=SbaseDop 
BSF 

doping uniform region=7 p.type conc=SbsfDop 








# Material properties 

Opaque contact [comment out for InGaP cell] 
material region=8 real.index=1.2 imag.index=1.8 

Vacuum (for zero reflection) [change to match window material (InGaP use 
Vacuum_AlInP) ] 





[for InGaP cell, comment out region 1] 
material region=1 index.file=Vacuum_InGaP.opt 
material region=2 index. file=VacuumIn20Ga80N. opt 
material region=3 index. file=VacuumIn20Ga80N. opt 

















#InGaN 
material material=InGaN EG300=2.6612 index.file=In20Ga80N.opt 





# Gold 
material material=Gold real.index=1.2 imag.index=1.8 








## Models [InGaP cell, 1; GaAs cell, 5&6; InGaNAs cell, 7] 
models region=1 CONMOB 
MODELS CHUANG CONMOB FLDMOB SRH OPTR PRINT 








## Light beams [GaAs b1,0.55-0.75,200 b2,0.75-0.88,65] 0.12-2.7,50 [630,825] 
beam num=1 x.origin=0 y.origin=S$lightY angle=90 back. refl 
power.file=AMOnrel.spec \ 

wavel.start=0.12 wavel.end=2.4 wavel.num=50 








struct outfile=SingleCell_webf.str 
#tonyplot SingleCell_webf.str 


solve init 
method gummel newton maxtraps=10 itlimit=25 
solve b1=0.9 


## Getting Isc for I-V curve points 
method newton maxtraps=10 itlimit=100 
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solve b1=0.95 

ract name="isc" 
isc=Sisc*$width3d 
il=Sisc/10 
12=$1i1+$isc/1 
13=$12+Sisc/1 
14=Si3+Sisc/1 
i15=$14+S$isc/1 
i16=Si5+Sisc/2 
i7=Si6+Sisc/2 
18=$i7+S$isc/2 
19=Si8+Sisc/2 
110=$i19+Sisc/20 
i111=$110+Sisc/40 
112=$111+$isc/40 
113=$112+S$isc/40 
114=$113+Sisc/40 
115=$114+Sisc/40 
116=$i15+Sisc/80 
i117=$i116+Sisc/80 
118=$117+S$isc/80 
119=$i18+Sisc/80 
120=$119+Sisc/80 


ex 
se 
se 
se 
se 
se 
se 
hs) 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
Ht 





CPO GE GE Ch.CE Gh cet cr Gh eh Cr Gheich “CF ck Geer “Gk. .cP. Cr CECE och cck ER oer. ct 


121=$1201 
122=$i214 
123=$1224 
124=$1234 
125=$1244 








0 
0 
0 
0 
0 
0 
0 
0 


max(i."cathode") 


tSisc/80 
tSisc/80 
tSisc/80 
t+Sisc/80 
tSisc/80-0.00001 


log outfile=In20Ga80N.log 


method newton maxtraps=10 itlimit=100 
solve b1=0.95 


contact name=anode current 
method newton maxtraps=10 itlimit=100 


Ht 


sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 





Pmax points 





























ve ianode=-$i25 bl= 
ve ianode=-$i24 bl= 
ve ianode=-$i23 bl= 
ve ianode=-$i22 bl= 
ve ianode=-$i21 bl= 
ve ianode=-$i20 bl= 
ve ianode=-$il9 bl= 
ve ianode=-$i18 bl= 
ve ianode=-$il17 bl= 
ve ianode=-$il6 bl= 
ve ianode=-$i15 bl= 
ve ianode=-$il14 bl= 
ve ianode=-$i13 bl= 
ve ianode=-$il2 bl= 
ve ianode=-$ill bl= 
ve ianode=-$i10 bl= 
ve ianode=-$i9 b1=0 
ve ianode=-$i8 b1=0 
ve ianode=-$i7 b1=0 
ve ianode=-$i6 b1=0 
ve ianode=-$i5 b1=0 
ve ianode=-$i4 b1=0 
ve ianode=-$i3 b1=0 





> 








Or COV: COON) CO. OC 'Or@rO.O:-¢ 


[InGaP 18-25; 


39D) 
~95 
95 
199) 
~95 
29S 
95 
399 
299 
95 
95 
cpa) 
95 
~95 
95 
299 


«95 
~95 
“93 
#95 
i995 
95 
~95 


GaAs 15-25; 
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InGaNAs 13-25; 


Ge 11-25] 


solve ianode=-$i2 b1=0.95 
solve ianode=-$il b1=0.95 











solve ianode=0 b1=0.95 

log off 

extract name="iv" curve(v."anode", i."cathode") outfile="IVcurvelIn20Gas0N.dat" 
tonyplot IVcurveIn20Ga80N.dat 


log outfile=doneIn20Ga80N.log 
log off 


B. SECOND JUNCTION: INo.57GAo.43N, EG=1.6 EV 


























go atlas 

set cellWidth=5.000000e+002 

set capWidthpercent=8.000000e+000 
set divs=1.000000e+001 

set contThick=1.000000e-001 

set capThick=3.000000e-001 

set capDop=1.000000e+020 

set windowThick=0.01 

set winDop=2.15e17 

set emitterThick=0.01 

#changed emitDop from lel6 to 1e20 
set emitDop=1lel16 

set baseThick=3.19467 

#changed basDop from 1lel16 to 1e20 
set baseDop=1lel16 

set bsfThick=0.03533 

set bsfDop=2.15e19 

set cellWidthDiv=S$cellWidth/Sdivs 
set width3d=100e6/$cellWidth 

set capWidth=0.01*ScapWidthpercent*S$cellWidth/2 
set capWidthDiv=ScapWidth/ (Sdivs/2) 
set cellWidthHalf=$cellWidth/2 

set bsfLo=0 

set bsfHi=SbsfLo-SbsfThick 

set bsfDiv=SbsfThick/Sdivs 

set baseLo=SbsfHi 

set baseHi=S$baseLo-SbaseThick 

set baseMid=SbaseLo-SbaseThick/2 
set baseDiv=SbaseThick/Sdivs 

set emitterLo=SbaseHi 

set emitterHi=SemitterLo-SemitterThick 
set emitterDiv=SemitterThick/S$divs 
set windowLo=SemitterHi 

set windowHi=SwindowLo-SwindowThick 
set windowDiv=SwindowThick/Sdivs 
set capLo=$windowHi 

set capHi=ScapLo-ScapThick 

#set capDiv=ScapThick/S$divs 

set contLo=S$capHi 

set contHi=$contLo-ScontThick 
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set contDiv=ScontThick/Sdivs 
set lightY=SemitterHi-5 


mesh width=Swidth3d 


## X-Mesh 

x.mesh loc=-$cellWidthHalf spac=$cellWidthDiv 

x.mesh loc=-ScapWidth spac=ScapWidthDiv 

x.mesh loc=ScapWidth spac=ScapWidthDiv 

x.mesh loc=ScellWidthHalf spac=$cellWidthDiv 
# Y—-Mesh 


Top contact 

-mesh loc=ScontHi spac=0 
y-.mesh loc=ScontLo spac=0 
Cap 
Window 

y-mesh loc=SwindowHi spac=SwindowDiv 
y-mesh loc=SwindowLo spac=SwindowDiv 
Emitter 

y.mesh loc=SemitterLo spac=SemitterDiv 
Base 
y.mesh loc=SbaseMid spac=SbaseDiv 
BSF 
y.mesh loc=SbsfHi spac=SbsfDiv 
y.mesh loc=SbsfLo spac=SbsfDiv 


ke 

















## Regions 

# Cap 

#region num=8 material=Vacuum x.min=-ScapWidth x.max=ScapWidth y.min=S$contHi 
y.max=ScontLo 

region num=1 material=InGaN x.min=-ScapWidth x.max=ScapWidth y.min=ScapHi 
y.max=ScapLo x.comp=0.57 

region num=2 material=Vacuum x.min=-S$cellWidthHalf x.max=—-ScapWidth 
y.min=ScontHi y.max=ScapLo 

region num=3 material=Vacuum x.min=ScapWidth x.max=ScellWidthHalf y.min=ScontHi 
y.max=ScapLo 

# Window [for Ge cell, use AlGaAs with x.comp=0.2] 






































region num=4 material=InGaN x.min=-ScellWidthHalf x.max=ScellWidthHalf 
y-.-min=SwindowHi y.max=SwindowLo x.comp=0.57 

region num=4 material=AlGaAs x.comp=0.2 x.min=-S$cellWidthHalf 
x.max=ScellWidthHalf y.min=SwindowHi y.max=SwindowLo 

Emitter 
region num=5 material=InGaN x.min=—-S$cellWidthHalf x.max=S$cellWidthHalf 
y-min=SemitterHi y.max=SemitterLo x.comp=0.57 

Base 
region num=6 material=InGaN x.min=—-S$cellWidthHalf x.max=ScellWidthHalf 
y.min=SbaseHi y.max=SbaseLo x.comp=0.57 

BSF 
region num=7 material=InGaN x.min=-ScellWidthHalf x.max=ScellWidthHalf 
y-.-min=SbsfHi y.max=SbsfLo x.comp=0.57 
## Electrodes [for InGaP cell, add cathode (gold) and remove cathode 
(conductor) ] 
electrode name=cathode material=Gold x.min=—-ScapWidth x.max=ScapWidth 
y.min=ScontHi y.max=ScontLo 
#electrode name=cathode x.min=-S$cellWidthHalf x.max=ScellWidthHalf 


y-.min=SwindowHi y.max=SwindowHi 
lectrod name=anod x.min=-ScellWidthHalf x.max=ScellWidthHalf y.min=S$bsfLo 
y.max=SbsfLo 
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Ht 


# Doping 


Cap 


doping uniform 


Window 





doping uniform 
Emitter 
doping uniform 


Base 


doping uniform 


BSF 


[for InGaP cell, 


region=1 n. 
region=4 n. 
region=5 n 
region=6 p. 
doping uniform region=7 p. 


# Material properties 
Opaque contact 


uncomment cap doping] 








Vacuum 


material region=2 index. 


[for 








GaAs 


InGaN 
material material=InGaN 


Gold 
material material=Gold real.index=1.2 imag.index=1.8 


odels 





material region=3 index. 








[InGaP cell, 1; 
models region=1 CONMOB 


MOD 
## Light beams 
beam num=1 


GaAs cell, 


x.origin=0 
power.file=PostJunctionl.spec.txt \ 
wavel.start=0.12 wavel.end=2.4 wavel.num=50 





struct outfile=SingleCell_webf.str 
#tonyplot SingleCell_webf.str 


solve init 
method gummel newton maxtraps=10 itlimit=25 
solve b1=0.9 


## Getting Isc for I-V curve points 


me 


solve b1=0.95 


ex 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 





ChGh ek (Gt GRoek “chaGh Cr:.ch. ct 


i5=$i44 
i6=$i5d4 
i7=$ié4 
i8=$i74 
i9=$i84 


110=$i9+Sisc/20 


tract name="isc" 


tSisc/10 
tSisc/20 
tSisc/20 
tSisc/20 
tSisc/20 





max(i."cathode") 
isc=Sisc*$width3d 
il=Sisc/10 

i12=Sil+Sisc/10 
i13=$i2+Sisc/10 
14=$i3+Sisc/10 


type conc=$capDop 

type conc=$winDop 
.type conc=SemitDop 
type conc=SbaseDop 
type conc=S$bsfDop 
[comment out for InGaP cell] 
material region=8 real.index=1.2 imag.index=1.8 


(for zero reflection) 
Vacuum_AlInP) ] 


[change to match window material (InGaP use 


InGaP cell, comment out region 1] 


file=VacuumiIn57Ga43N.opt 
file=VacuumiIn57Ga43N.opt 


5&6; InGaNAs cell, 


ELS CHUANG CONMOB FLDMOB SRH OPTR PRINT 


[GaAs b1,0.55-0.75,200 b2,0.75-0.88, 65] 


y.origin=SlightyY 


thod newton maxtraps=10 itlimit=100 
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EG300=1.559 index. file=In57Ga43N. opt 


7] 


O.12=2..7750 
angle=90 


[630,825] 
back.refl 


se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
se 
Ht 





CF Gk Gh cr Gk Ach. (Ch CECE 2G “Er SER GE Gh oct 


i111=$110+Sisc/40 
112=$111+$isc/40 
113=$112+S$isc/40 
114=$i113+Sisc/40 
115=$114+Sisc/40 
116=$115+Sisc/80 
117=$1i16+Sisc/80 
118=$117+Sisc/80 
119=$i18+S$isc/80 
120=$i119+Sisc/80 


121=$1204 
122=$1214 
123=$1224 
124=$1234 
125=$1244 








t+Sisc/80 
t+Sisc/80 
tSisc/80 
t+Sisc/80 
tSisc/80-0.00001 


log outfile=In57Ga43N.log 


method newton maxtraps=10 itlimit=100 
solve b1=0.95 


contact name=anode current 
method newton maxtraps=10 itlimit=100 












































## Pmax points 

solve ianode=-$i25 bl= 
solve ianode=-$i24 bl= 
solve ianode=-$i23 bl= 
solve ianode=-$i22 bl= 
solve ianode=-$i21 bl= 
solve ianode=-$i20 bl= 
solve ianode=-$il19 bl= 
solve ianode=-$i18 bl= 
solve ianode=-$il7 bl= 
solve ianode=-$il6 bl= 
solve ianode=-$il15 bl= 
solve ianode=-$il14 bl= 
solve ianode=-$il3 bl= 
solve ianode=-$il2 bl= 
solve ianode=-$ill bl= 
solve ianode=-$i10 bl= 
solve ianode=-$i9 b1=0 
solve ianode=-$i8 b1=0 
solve ianode=-$i7 b1=0 
solve ianode=-$i6 b1=0 
solve ianode=-$i5 b1=0 
solve ianode=-$i4 b1=0 
solve ianode=-$i3 b1=0 
solve ianode=-$i2 b1=0 
solve ianode=-$il b1=0 
## 

solve ianode=0 b1=0.95 
log off 


extract name="iv" 


tonyplot IVcurveIn57Ga43N.dat 


log outfile=doneIn57Ga43N.log 


log off 


> 








OC O°@ OO O'O"O::C-O'O'O Ore O:-°* 


[InGaP 18-25; 


295 
~95 
~95 
95 
+99 
95 
95 
~95 
295 
~95 
~95 
~95 
95 
ac) 
95 
95 


+95 
95 
Paice) 
95 
“9D 
95 
95 
~95 
«95 


GaAs 15-25; InGaNAs 13-25; Ge 11-25] 


curve (v."anode", i."cathode") 
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outfile="IVcurvelIn57Ga43N.dat" 


Cc. THIRD JUNCTION: INo.¢6sGAo.32N, EG=1.31 EV 


go atlas 


set cellWidth=5.000000e+002 

set capWidthpercent=8.000000e+000 
set divs=1.000000e+001 

set contThick=1.000000e-001 


set capThick=3.000000e-001 


set capDop=1.000000e+020 
set windowThick=0.01 

set winDop=2.15e17 

s mitterThick=0.01 





#changed emitDop from lel6 to 1e20 
set emitDop=1lel16 

set baseThick=3.19467 

#changed basDop from lel6é to 1e20 























set baseDop=1el16 

set bsfThick=0.03533 

set bsfDop=2.15e19 

set cellWidthDiv=S$cellWidth/Sdivs 
set width3d=100e6/$cellWidth 

set capWidth=0.01*ScapWidthpercent*S$cellWidth/2 
set capWidthDiv=ScapWidth/ (Sdivs/2) 
set cellWidthHalf=$cellWidth/2 

set bsfLo=0 

set bsfHi=SbsfLo-SbsfThick 

set bsfDiv=SbsfThick/Sdivs 

set baseLo=SbsfHi 

set baseHi=S$baseLo-SbaseThick 

set baseMid=SbaseLo-SbaseThick/2 
set baseDiv=SbaseThick/Sdivs 

set emitterLo=SbaseHi 

set emitterHi=SemitterLo-SemitterThick 
set emitterDiv=SemitterThick/S$divs 
set windowLo=SemitterHi 

set windowHi=SwindowLo-SwindowThick 
set windowDiv=SwindowThick/Sdivs 
set capLo=$windowHi 

set capHi=ScapLo-ScapThick 


#set capDiv=ScapThick/S$divs 


t contLo=ScapHi 
set contHi=$contLo-ScontThick 
t contDiv=$contThick/S$divs 








set lightY=SemitterHi-5 


mesh width=$width3d 

## X-Mesh 

x.mesh loc=-$cellWidthHalf spac=$cellWidthDiv 
x.mesh loc=—-$capWidth spac=ScapWidthDiv 
x.mesh loc=ScapWidth spac=ScapWidthDiv 

x.mesh loc=ScellWidthHalf spac=$cellWidthDiv 
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# Y—-Mesh 

Top contact 

y-mesh loc=ScontHi spac=0 
y-mesh loc=ScontLo spac=0 
Cap 
Window 

y-mesh loc=SwindowHi spac=SwindowDiv 
y-mesh loc=SwindowLo spac=SwindowDiv 
Emitter 

y-mesh loc=SemitterLo spac=SemitterDiv 
Base 
y.mesh loc=SbaseMid spac=SbaseDiv 














y-mesh loc=SbsfHi spac=SbsfDiv 
y.mesh loc=SbsfLo spac=SbsfDiv 





## Regions 

# Cap 

#region num=8 material=Vacuum x.min=—-ScapWidth x.max=ScapWidth y.min=ScontHi 
y.max=ScontLo 

region num=1 material=InGaN x.min=-ScapWidth x.max=ScapWidth y.min=ScapHi 
y.max=ScapLo x.comp=0.68 

region num=2 material=Vacuum x.min=-S$cellWidthHalf x.max=—-ScapWidth 
y.min=ScontHi y.max=ScapLo 

region num=3 material=Vacuum x.min=ScapWidth x.max=ScellWidthHalf y.min=ScontHi 
y.max=ScapLo 

# Window [for Ge cell, use AlGaAs with x.comp=0.2] 



































region num=4 material=InGaN x.min=—-S$cellWidthHalf x.max=S$cellWidthHalf 
y-.-min=SwindowHi y.max=SwindowLo x.comp=0.68 

region num=4 material=AlGaAs x.comp=0.2 x.min=-S$cellWidthHalf 
x.max=ScellWidthHalf y.min=SwindowHi y.max=SwindowLo 

Emitter 
region num=5 material=InGaN x.min=-ScellWidthHalf x.max=ScellWidthHalf 
y-min=SemitterHi y.max=SemitterLo x.comp=0.68 

Base 
region num=6 material=InGaN x.min=-ScellWidthHalf x.max=ScellWidthHalf 
y.min=SbaseHi y.max=SbaseLo x.comp=0.68 

BSF 
region num=7 material=InGaN x.min=-ScellWidthHalf x.max=ScellWidthHalf 
y-.min=SbsfHi y.max=SbsfLo x.comp=0.68 
## Electrodes [for InGaP cell, add cathode (gold) and remove cathode 
(conductor) ] 
electrode name=cathode material=Gold x.min=-ScapWidth x.max=ScapWidth 
y.min=ScontHi y.max=ScontLo 
#electrode name=cathode x.min=—-ScellWidthHalf x.max=ScellWidthHalf 


y.min=SwindowHi y.max=SwindowHi 
lectrod name=anod x.min=-ScellWidthHalf x.max=$cellWidthHalf y.min=SbsfLo 
y.max=SbsfLo 








# Doping [for InGaP cell, uncomment cap doping] 
Cap 

doping uniform region=1 n.type conc=ScapDop 
Window 

doping uniform region=4 n.type conc=$winDop 
Emitter 

doping uniform region=5 n.type conc=SemitDop 
Base 








doping uniform region=6 p.type conc=SbaseDop 
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# BSF 
doping uniform region=7 p.type conc=SbsfDop 


## Material properties 

# Opaque contact [comment out for InGaP cell] 
#material region=8 real.index=1.2 imag.index=1.8 
# Vacuum (for zero reflection) 

material region=2 index. file=VacuumIn68Ga32N.opt 
material region=3 index. file=VacuumIn68Ga32N.opt 











#InGaN 
material material=InGaN EG300=1.3068 index. file=In68Ga32N.opt 





# Gold 
material material=Gold real.index=1.2 imag.index=1.8 








## Models [InGaP cell, 1; GaAs cell, 5&6; InGaNAs cell, 7] 
models region=1 CONMOB 
MODELS CONMOB FLDMOB SRH OPTR PRINT 





## Light beams [GaAs b1,0.55-0.75,200 b2,0.75-0.88,65] 0.12-2.7,50 [630,825] 
beam num=1 x.origin=0 y.origin=S$lightY angle=90 back. refl 
power.file=PostJunction2.spec.txt \ 

wavel.start=0.12 wavel.end=2.4 wavel.num=50 


struct outfile=SingleCell_webf.str 
#tonyplot SingleCell_webf.str 





solve init 
method gummel newton maxtraps=10 itlimit=25 
solve b1=0.9 


## Getting Isc for I-V curve points 











method newton maxtraps=10 itlimit=100 
solve b1=0.95 

extract name="isc" max(i."cathode") 
set isc=Sisc*$width3d 
set il=Sisc/10 

set i2=$11+Sisc/10 
set i3=$12+Sisc/10 
set i4=$13+Sisc/10 
set i5=$14+Sisc/10 
set i6=$i5+Sisc/20 
set i7=$16+Sisc/20 
set i8=$1i17+Sisc/20 
set 19=$i8+Sisc/20 
set 110=$19+Sisc/20 
set i11=$110+Sisc/40 
set i12=$111+Sisc/40 
set 113=$112+Sisc/40 
set i14=$113+Sisc/40 
set 115=$i14+Sisc/40 
set i116=$i15+$isc/80 
set i117=$i16+Sisc/80 
set i118=$i17+Sisc/80 
set i119=$i18+$isc/80 
set i120=$i19+$isc/80 
set i21=$120+Sisc/80 
set i22=$121+Sisc/80 
set i23=$122+Sisc/80 
set i24=$123+Sisc/80 


719 


set 125=$124+Sisc/80-0.00001 
Ht 


log outfile=In68Ga32N.log 


method newton maxtraps=10 itlimit=100 
solve b1=0.95 


contact name=anode current 
method newton maxtraps=10 itlimit=100 


## Pmax points [InGaP 18-25; GaAs 15-25; InGaNAs 13-25; Ge 11-25] 















































solve ianode=-$i25 b1=0.95 
solve ianode=-$i24 b1=0.95 
solve ianode=-$i23 b1=0.95 
solve ianode=-$i22 b1=0.95 
solve ianode=-$i21 b1=0.95 
solve ianode=-$i20 b1=0.95 
solve ianode=-$il9 b1=0.95 
solve ianode=-$il8 b1=0.95 
solve ianode=-$i1l7 b1=0.95 
solve ianode=-$il6 b1=0.95 
solve ianode=-$il5 b1=0.95 
solve ianode=-$il14 b1=0.95 
solve ianode=-$il3 b1=0.95 
solve ianode=-$il2 b1=0.95 
solve ianode=-$ill b1=0.95 
solve ianode=-$il0 b1=0.95 
solve ianode=-$i9 b1=0.95 
solve ianode=-$i8 b1=0.95 
solve ianode=-$i7 b1=0.95 
solve ianode=-$i6 b1=0.95 
solve ianode=-$i5 b1=0.95 
solve ianode=-$i4 b1=0.95 
solve ianode=-$i3 b1=0.95 
solve ianode=-$i2 b1=0.95 
solve ianode=-$il b1=0.95 
## 

solve ianode=0 b1=0.95 





log off 
extract name="iv" curve(v."anode", i."cathode") outfile="IVcurvelIn68Ga32N.dat" 
tonyplot IVcurveIn68Ga32N.dat 


log outfile=doneIn68Ga32N.log 
log off 


D. BOTTOM JUNCTION: INo.7sGAo.22N, EG=1.11 EV 


go atlas 


set cellWidth=5.000000e+002 

set capWidthpercent=8.000000e+000 
set divs=1.000000e+001 

set contThick=1.000000e-001 

set capThick=3.000000e-001 

set capDop=1.000000e+020 
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set windowThick=0.01 

set winDop=2.15e17 

set emitterThick=0.01 

#changed emitDop from 1lel6 to 1e20 
set emitDop=lel6 

set baseThick=3.19467 

#changed basDop from 1lel16 to 1e20 
set baseDop=1el16 

set bsfThick=0.03533 

set bsfDop=2.15e19 

set cellWidthDiv=S$cellWidth/Sdivs 
set width3d=100e6/$cellWidth 

set capWidth=0.01*S$capWidthpercent*S$cellWidth/2 
set capWidthDiv=ScapWidth/ (Sdivs/2) 
set cellWidthHalf=$cellWidth/2 

set bsfLo=0 

set bsfHi=SbsfLo-SbsfThick 

set bsfDiv=SbsfThick/Sdivs 

set baseLo=SbsfHi 

set baseHi=S$baseLo-SbaseThick 

set baseMid=SbaseLo-SbaseThick/2 
set baseDiv=SbaseThick/Sdivs 

set emitterLo=SbaseHi 

set emitterHi=SemitterLo-SemitterThick 
set emitterDiv=SemitterThick/S$divs 
set windowLo=SemitterHi 

set windowHi=SwindowLo-SwindowThick 
set windowDiv=SwindowThick/Sdivs 
set capLo=$windowHi 

set capHi=ScapLo-ScapThick 

#set capDiv=ScapThick/S$divs 

set contLo=$capHi 

set contHi=$contLo-S$contThick 

set contDiv=ScontThick/Sdivs 

set lightY=SemitterHi-5 





mesh width=$width3d 

## X-Mesh 

x.mesh loc=-$cellWidthHalf spac=$cellWidthDiv 
x.mesh loc=—-$capWidth spac=ScapWidthDiv 
x.mesh loc=ScapWidth spac=ScapWidthDiv 

x.mesh loc=ScellWidthHalf spac=$cellWidthDiv 


# Y—-Mesh 

Top contact 

-mesh loc=ScontHi spac=0 
y-mesh loc=ScontLo spac=0 
Cap 
Window 

y-mesh loc=SwindowHi spac=SwindowDiv 
y-mesh loc=SwindowLo spac=SwindowDiv 
Emitter 

y.mesh loc=SemitterLo spac=SemitterDiv 
Base 


kK 

















81 


-mesh loc=SbaseMid spac=S$baseDiv 
BSF 

-mesh loc=SbsfHi spac=S$bsfDiv 
-mesh loc=SbsfLo spac=S$bsfDiv 


KEK SEA 


## Regions 

# Cap 

#region num=8 material=Vacuum x.min=—-ScapWidth x.max=ScapWidth y.min=ScontHi 
y.max=ScontLo 

region num=1 material=InGaN x.min=-ScapWidth x.max=ScapWidth y.min=ScapHi 
y.max=ScapLo x.comp=0.78 

region num=2 material=Vacuum x.min=-S$cellWidthHalf x.max=—-ScapWidth 
y.min=ScontHi y.max=ScapLo 

region num=3 material=Vacuum x.min=ScapWidth x.max=ScellWidthHalf y.min=ScontHi 
y.max=ScapLo 

# Window [for Ge cell, use AlGaAs with x.comp=0.2] 



































region num=4 material=InGaN x.min=—-S$cellWidthHalf x.max=ScellWidthHalf 
y-.min=SwindowHi y.max=SwindowLo x.comp=0.78 

region num=4 material=AlGaAs x.comp=0.2 x.min=-S$cellWidthHalf 
x.max=ScellWidthHalf y.min=SwindowHi y.max=SwindowLo 

Emitter 
region num=5 material=InGaN x.min=-ScellWidthHalf x.max=ScellWidthHalf 
y-min=SemitterHi y.max=SemitterLo x.comp=0.78 

Base 
region num=6 material=InGaN x.min=-ScellWidthHalf x.max=ScellWidthHalf 
y.min=SbaseHi y.max=SbaseLo x.comp=0.78 

BSF 
region num=7 material=InGaN x.min=-ScellWidthHalf x.max=ScellWidthHalf 
y-.min=SbsfHi y.max=SbsfLo x.comp=0.78 
## Electrodes [for InGaP cell, add cathode (gold) and remove cathode 
(conductor) ] 
electrode name=cathode material=Gold x.min=—-ScapWidth x.max=ScapWidth 
y.min=ScontHi y.max=ScontLo 
#electrode name=cathode x.min=—-S$cellWidthHalf x.max=ScellWidthHalf 


y.min=SwindowHi y.max=SwindowHi 
lectrod name=anod x.min=-ScellWidthHalf x.max=$cellWidthHalf y.min=SbsfLo 
y.max=SbsfLo 





# Doping [for InGaP cell, uncomment cap doping] 











Cap 

doping uniform region=1 n.type conc=ScapDop 
Window 

doping uniform region=4 n.type conc=$winDop 
Emitter 

doping uniform region=5 n.type conc=SemitDop 
Base 

doping uniform region=6 p.type conc=SbaseDop 
BSF 

doping uniform region=7 p.type conc=SbsfDop 


# Material properties 

Opaque contact [comment out for InGaP cell] 
material region=8 real.index=1.2 imag.index=1.8 
Vacuum (for zero reflection) 
material region=2 index.file=VacuumInGa22N.opt 
material region=3 index. file=VacuumInGa22N.opt 

















#InGaN 
material material=InGaN EG300=1.1076 index.file=InGa22N.opt 
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# Gold 
material material=Gold real.index=1.2 imag.index=1.8 





## Models [InGaP cell, 1; GaAs cell, 5&6; InGaNAs cell, 7] 
models region=1 CONMOB 
MODELS CONMOB FLDMOB SRH OPTR PRINT 





## Light beams [GaAs b1,0.55-0.75,200 b2,0.75-0.88,65] 0.12-2.7,50 [630,825] 
beam num=1 x.origin=0 y.origin=SlightY angle=90 back.refl 
power.file=PostJunction3.spec.txt \ 

wavel.start=0.12 wavel.end=2.4 wavel.num=50 





struct outfile=SingleCell_webf.str 
#tonyplot SingleCell_webf.str 


solve init 
method gummel newton maxtraps=10 itlimit=25 


solve b1=0.9 


## Getting Isc for I-V curve points 











method newton maxtraps=10 itlimit=100 
solve b1=0.95 

extract name="isc" max(i."cathode") 
set isc=Sisc*$width3d 

set il=Sisc/10 

set i2=$1i11+Sisc/10 

set i3=$12+Sisc/10 

set i4=$13+Sisc/10 

set i5=$14+Sisc/10 

set i6=$i5+Sisc/20 

set i7=$16+Sisc/20 

set i8=$1i7+Sisc/20 

set i19=$i8+Sisc/20 

set 110=$19+Sisc/20 

set i11=$110+Sisc/40 

set i12=$111+Sisc/40 

set 113=$112+Sisc/40 

set i14=$113+Sisc/40 

set 115=$i14+Sisc/40 

set i116=$i115+$isc/80 

set i117=$i16+Sisc/80 

set i118=$i17+Sisc/80 

set i119=$i18+$isc/80 

set i120=$i19+$isc/80 

set i21=$120+Sisc/80 

set i22=$121+Sisc/80 

set i23=$122+Sisc/80 

set i24=$123+Sisc/80 

set 125=$124+$isc/80-0.00001 
## 


log outfile=In78Ga22N.log 


method newton maxtraps=10 itlimit=100 
solve b1=0.95 


contact name=anode current 
method newton maxtraps=10 itlimit=100 


## Pmax points [InGaP 18-25; GaAs 15-25; InGaNAs 13-25; Ge 11-25] 
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solve ianode=-$i25 b1=0.95 
solve ianode=-$i24 b1=0.95 
solve ianode=-$i23 b1=0.95 
solve ianode=-$i22 b1=0.95 
solve ianode=-$i21 b1=0.95 
solve ianode=-$i20 b1=0.95 
solve ianode=-$il9 b1=0.95 
solve ianode=-$il8 b1=0.95 
solve ianode=-$i1l7 b1=0.95 
solve ianode=-$il6 b1=0.95 
solve ianode=-$il5 b1=0.95 
solve ianode=-$il4 b1=0.95 
solve ianode=-$il3 b1=0.95 
solve ianode=-$il2 b1=0.95 
solve ianode=-$ill b1=0.95 
solve ianode=-$il0 b1=0.95 
solve ianode=-$i9 b1=0.95 
solve ianode=-$i8 b1=0.95 
solve ianode=-$i7 b1=0.95 
solve ianode=-$i6 b1=0.95 
solve ianode=-$i5 b1=0.95 
solve ianode=-$i4 b1=0.95 
solve ianode=-$i3 b1=0.95 
solve ianode=-$i2 b1=0.95 
solve ianode=-$il b1=0.95 
## 

solve ianode=0 b1=0.95 





log off 
extract name="iv" curve(v."anode", i."cathode") outfile="IVcurvelIn78Ga22N.dat" 
tonyplot IVcurveIn78Ga22N.dat 


log outfile=doneIn78Ga22N.log 
log off 
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APPENDIX B: MATLAB CODE 


The Matlab code provided in this appendix provided 


auxiliary support in the interpretation and display Sivaco 





Atlas log files. 





Original code for some of the Matlab functions and 


scripts was obtained from Michalopolous [1], Bates, [2], 





Green [3], and Canfield [4]. Modifications to the code were 


made to support this thesis. 


A. INDIUM GALLIUM NITRIDE BAND GAP CALCULATIONS 


SInGaN Bandgap calculations 
SBaldomero Garcia 
GLa: 


SConstants 


EgGaN=3.42; 
EgInN=0.77; 

b=1.43; 
x=[0:0.0001:1]; 
x=.67; 
InConcentration=x 
GaConcentration=1-x 





ae ole 


oe 





SFormula for finding EgInGaN 
EgInGaN=((1-x) .*EgGaN) + (x.*EgInN) —-(b.*x.* (1-x) ) 














sFormula for finding x 
solve (' (( (1-x) *EgGaN) + (x*EgInN) -— (b*x* (1-x) )) -EgInGaN', 'x"'); 

















% The other EgInGaN that we're solving for is EgInGaN=1.1 
% EgInGaN=2.095 

% Gacomp=1-(1/2/b* (EgGaN-EgInN+b- (EgGaN*2-2*EgInN*EgGaN- 
2*b*EgGaNtEgInN*2-2*EgInN*b+b*2+4*b*EgInGaN) *(1/2))) 

% Incomp=(1/2/b* (EgGaN-EgInN+b- (EgGaN*2-2*EgInN*EgGaN- 
2*b*EgGaNtEgInN*2-2*EgInN*b+b*2+4*b*EgInGaN) *(1/2))) 






















































































S$Plot 

plot (x,EgInGaN, 'b'); 

xlabel('Indium concentration (unitless ratio)'); 
ylabel('InGaN Band Gap (eV)'); 











title('Indium concentration vs InGaN Band Gap'); 
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B. CONVERSION FROM DIELECTRIC CONSTANTS (EPSILONS) TO 
REFRACTION COEFFICIENTS (N, K) 


Optical data needed to be provided n-k format. Some of 


the data obtained was available in epsilon form. Therefore, 





a conversion needed to be made. Michalopoulos generated the 


function below [1]. 


oe 





E2NK Convert th 1, @2 pairs into n, k pairs’. 
& [n k] = E2NK(el, e2) 








% (c)2001 by P. Michalopoulos 








function [n,k] = e2nk(el, e2) 
ap = (el + sqrt(el.*2 + e2.%2)) / 2; 
an = (el - sqrt(el.*2 + e2.%2)) / 2; 
app = ap >= 0; 
anp = an >= 0; 
err = (app < 0) & (anp < 0); 
err = sum(err); 
if err ~= 0 
disp ('ERROR!') 
end 
a= ap .* app + an .* anp; 
n = sqrt(a); 
k = e2 ./ (2 * n); 
Cc. CONVERSION FROM PHOTON ENERGY (EV) TO WAVELENGTH (UM) 





Data was found in photon energy and wavelength form. 





In order to make comparisons, conversion from one form to 





another was necessary. Functions by Michalopoulos [1] are 


given below. 





% EV2UM Converts photon energy (eV) into wavelength (um). 


% (c)2001 by P. Michalopoulos 


function um = ev2um(ev) 
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h = 6.6260755e-34; 
eV = 1.60218e-19; 
c = 2.99792458e8; 
ev = ev * eV; 

f= eyv-/ hy 

wavel =c ./ f; 

um = wavel / le-6; 


oe 


UM2EV 





oe 


(CG) 2001 by PB. 


function ev = 


Michalopoulos 


um2ev (um) 


Converts photon wavelength 


(um) 


into 


energy (eV). 


h = 6.6260755e-34; 

eV = 1.60218e-19; 

c = 2.99792458e8; 

wavel = um * le-6; 

f =c ./ wavel; 

ev =h * f; 

ev = ev ./ eV; 

D. IV CURVE PLOTS FOR INDIUM GALLIUM NITRIDE QUAD 


JUNCTION SOLAR CELL 


% Code modified from Canfield's source code 

Data obtained from log files from InGaN quad-junction solar cell 
% Baldomero Garcia 
ele: 
6 In20Ga80N with AMO spectrum 
current1l=[0.0176131398 0.0176031407 





0.0175505986 0.0174980573 


0.0173929742 0.0171728094 0.0169526469 0.0167324843 
0.0165123158 0.0162921532 0.0160719881 0.0158518253 
0.0156316602 0.0154114983 0.0149711686 0.0145308402 
0.0140905121 0.0136501825 0.0132098541 0.0123291993 
0.0114485415 0.0105678837 0.0096872266 0.0088065691 
0.0070452558 0.0052839419 0.0035226276 0.0017613142 
0.0000000001]; 

voltagel=[0.0000000000 0.2970913579 1.1822462140 1.5075176680 
1.7472436590 1.8900757790 1.9378145180 1.9663764170 
1.9875670060 2.0047710060 2.0194136280 2.0322369610 
2.0436835120 2.0540360290 2.0722008890 2.0877645130 
2.1013319790 2.1132957310 2.1239292560 2.1419700890 
2.1566166780 2.1686629450 2.1786938920 2.1871519430 
2.2006132650 2.2108811030 2.2190289820 2.2257052930 
2.2313181110]; 


oe 


In57Ga43N with postjunctionl spectrum 
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current2=[0.0171128720 


0.0171028783 





0.0168989671 


0.0166850573 


0.0164711457 .0162572338 0.0160433258 0.0158294126 
0.0156155026 -0154015923 0.0151876797 0.0149737698 
0.0145459485 .0141181240 0.0136903032 0.0132624807 
0.0128346588 .0119790166 0.0111233713 0.0102677280 
0.0094120828 -0085564398 0.0068451513 0.0051338641 
0.0034225761 .0017112879 -0.0000000001]; 
voltage2=[0.0000000000 0.1305081858 0.8192581792 0.9015131674 
0.9331338594 - 9534997163 0.9689573887 0.9815965012 
0.9923675261 .0017892410 1.0101786220 1.0177449700 
1.0309568480 - 0422028240 1.0519502020 1.0605082500 
1.0680958050 .0809744750 1.0915085520 1.1002972060 
1.1077533400 .1141715070 1.1247003520 1.1330384420 
1.1398643710 . 1456000250 .1505207070]; 

In68Ga32N with postjunction2 spectrum 
current3=[0.0141231105 0.0141131184 0.0140298496 0.0139465807 
0.0137700408 .0135935022 0.0134169620 0.0132404230 
0.0130638854 -0128873465 0.0127108077 0.0125342689 
0.0123577303 .0120046521 0.0116515739 0.0112984952 
0.0109454185 .0105923395 0.0098861822 0.0091800273 
0.0084738721 .0077677158 0.0070615601 0.0056492474 
0.0042369367 -0028246240 0.0014123119 0.0000000001]; 
voltage3=[0.0000000000 0.1311096030 0.5441673681 0.6320277025 
0.6842499725 . 7083967768 0.7248954021 0.7377135412 
0.7483143199 0.7574074704 0.7653951663 0.7725296065 
0.7789805527 0.7902824668 0.7999460470 0.8083630634 
0.8157932979 0.8224203640 0.8337763994 0.8431914560 
0.8511524069 0.8579918718 0.8639469242 0.8738573284 
0.8818299562 0.8884345077 0.8940340114 0.8988707325]; 
6 In78Ga22N with potjunction3 spectrum 














current4=[0.0128841840 0.0128741794 0.0127231275 0.0125620760 
0.0124010228 0122399704 0.0120789198 0.0119178658 
0.0117568129 0115957622 0.0114347108 0.0112736577 
0.0109515531 0106294488 0.0103073439 0.0099852400 
0.0096631342 0090189240 0.0083747191 0.0077305090 
0.0070862990 .0064420895 0.0051536706 0.0038652545 
0.0025768364 0012884183 -0.0000000001]; 
voltage4=[0.0000000000 0.1078032006 0.4688474643 0.5100668879 
0.5309258388 - 5454971566 0.5569116230 0.5663876299 
0.5745324881 0.5816962863 0.5881009992 0.5938970077 
0.6040640671 0.6127739225 0.6203776540 0.6271079482 
0.6331288820 0.6434973415 0.6521554853 0.6595301586 
0.6659109540 0.6715037308 0.6808929507 0.6885228391 
0.6948948317 0.7003322584 0.7050534598]; 

% Combined 

current5=[0 0.001991 0.004011 0.006003 0.007994 0.01001 
0.01178 0.01198 0.0126 0.01275 0.01286 0.0128841840]; 
voltage5=[4.9793 4.9542 4.9195 4.8744 4.8185 4.7271 4.5995 
4.5723 4.4736 4.3906 4.0906 0]; 


Pl=currentl.*voltagel; 
[a b]=max(P1); 
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Pmaxl=voltagel (b) *currentl (b); 
Vocl=max (voltagel) ; 

Iscl=max (currentl1)j; 

FF1=Pmaxl1/ (Vocl1*Iscl) 
Ef£1=100*Pmaxl1/ (.1353) 


P2=current2.*voltage2; 

[c d]=max(P2); 
Pmax2=voltage2 (d) *current2 (d) ; 
Voc2=max (voltage2) ; 

Isc2=max (current2); 

FF2=Pmax2/ (Voc2*Isc2) 
Eff£2=100*Pmax2/ (.1353) 





P3=current3.*voltage3; 

[e f]=max(P3); 
Pmax3=voltage3 (f) *current3 (f); 
Voc3=max (voltage3) ; 

Isc3=max (current3); 

FF3=Pmax3/ (Voc3*Isc3) 
Ef£3=100*Pmax3/ (.1353) 





P4=current4.*voltage4; 

[g h]=max(P4); 
Pmax4=voltage4 (h) *currenté4 (h); 
Voc4=max (voltage4) ; 

Isc4=max (current4) ; 

FF4=Pmax4/ (Voc4*Isc4) 
Ef£f£4=100*Pmax4/ (.1353) 





P5=current5.*voltage5; 
[k l]=max(P5); 
Pmax5=voltage5(1)*current5 (1); 
Voc5=max (voltageS5) ; 

Isc5=max (currentd5); 

FF5=Pmax5/ (Voc5*Isc5) 
Ef£5=100*Pmax5/ (.1353) 











plot (voltagel,currentl,'r'); 
sgrid on; 

hold on; 

splot (voltagel (b),currentl(b), 
hold on; 

plot (voltage2,current2,'r'); 
hold on; 

splot (voltage2 (d),current2(d), 
hold on; 

plot (voltage3,current3,'r'); 
hold on; 

splot (voltage3(f),current3(f), 
hold on; 

plot (voltage4,current4, 'r'); 
hold on; 

splot (voltage4 (h),current4(h), 
hold on; 








"o'); 
NOT ie 
"o'); 
ok 


lot (voltage5,current5,'b'); 
old on; 
ot (voltage5(1),current5(1),'or'); 
old off; 

label ('Voltage (V)') 
ylabel('Current Density (A/cm%*2)') 
title('Solar Cell IV Curve') 
axis([0 5.5 0 0.020]) 














x 7'O D'O 

















E. AIR MASS ZERO PLOTS 


Air Mass Zero (AMO) data was obtained from NREL [19] 


and plotted using the Matlab function provided by 








Michalopoulos [1]. It was modified to provide an additional 


plot. 
% AMSPECTRUMS Plot the data contained in spectrums.mat 





%& Check the source code and uncomment the required lines. 

% (c)2001 by P. Michalopoulos 

h = 6.6260755e-34; 

1.60218e-19; 

c = 2.99792458e8; 

load('Data\spectrums') 

plot (AMO_wavel, AMO_intl, 'b', AM15a_wavel, AM15a_intl, '‘r'), grid on 
title('AM 0 and AM 1.5 spectra'), xlabel('Wavelength [um]'), 
ylabel('Irradiance [W / cm*2*um]') 

axis([0 3 0 0.23]) 





0) 
s 
ll 








The code above was modified to display AMO with respect 
to photon energy and wavelength, respectively. 


Figure (1) 
plot (um, iradiance,'r'); 


























xlabel('Wavelength (um) ') 
ylabel ('Irradiance W/cm*2*um') 
title ('AMO') 

axis([0 2.5 0 0.25]) 

Figure (2) 

plot (ev,iradiance,'b'); 
xlabel('Energy (eV) ') 

ylabel ('Irradiance W/cm*2um') 
title ('AMO') 

axis([0 6.0 0 0.25]) 
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